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v  ABSTRACT 

This  report  considers  plane-vave  scattering  by  a  homogeneous 
material  body  of  revolution.  The  problem  is  formulated  in  terms  of 
equivalent  electric  and  magnetic  currents  over  the  surface  which 
defines  the  body.  Application  of  the  boundary  conditions  leads  to 
four  simultaneous  surface  integral  equations  to  be  satisfied  by  the 
two  unknown  equivalent  currents,  electric  and  magnetic.  The  set  of 
four  equations  is  reduced  to  a  coupled  pair  of  equations  by  taking 
linear  combinations  of  the  original  four  equations.  Because  many 
pairs  of  linear  combinations  are  possible,  there  are  many  surface 
integral  equation  formulations  for  the  problem.  Two  formulations 
commonly  encountered  in  the  literature  are  discussed  and  solved  by 
the  method  of  moments.  A  general  computer  program  for  material 
bodies  of  revolution  is  developed,  listed,  and  documented.  Examples 
of  numerical  computations  are  given  for  dielectric  spheres  and  a 
finite  dielectric  cylinder.  The  computed  results  for  the  sphere  are 
compared  to  the  exact  series  solution  obtained  by  separation  of 
variables , 
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PART  ONE 


ELECTROMAGNETIC  SCATTERING  FROM  A  HOMOGENEOUS  MATERIAL 
BODY  OF  REVOLUTION 

THEORY  AND  EXAMPLES 


I .  INTRODUCTION 

The  problem  of  plane-wave  scattering  by  a  homogeneous  material  body 
of  revolution  is  formulated  in  terms  of  equivalent  electric  and  magnetic 
currents  over  the  body  surface.  Application  of  boundary  conditions  leads 
to  a  set  of  four  integral  equations  to  be  satisfied.  Linear  combinations 
of  these  four  equations  lead  to  a  coupled  pair  of  equations  to  be  solved. 
One  choice  of  combination  constants  gives  the  formulation  described  by 
Poggio  and  Miller  [1].  This  formulation  has  been  applied  to  material 
cylinders  by  Chang  and  Harrington  [2],  and  to  material  bodies  of  revolu¬ 
tion  by  Wu  [3].  We  will  call  this  choice  the  PMCHW  formulation  (formed 
by  the  initials  of  the  above  cited  investigators). 

Another  choice  of  combination  constants  gives  the  formulation 
obtained  by  Miiller  [4].  This  formulation  has  been  applied  to  dielectric 
cylinders  by  Solodukhov  and  Vasil'ev  [5]  and  by  Morita  [6],  and  to  bodies 


[1]  A.  J.  Poggio  and  E.  K.  Miller,  "Integral  Equation  Solutions  of  Three- 
dimensional  Scattering  Problems,"  Chap.  4  of  Computer  Techniques  for 
Electromagnetics ,  edited  by  R.  Mittra,  Pergamon  Press,  1973,  Equa¬ 
tion  (4.17). 

[2]  Yu  Chang  and  R.  F.  Harrington,  "A  Surface  Formulation  for  Characteristic 
Modes  of  Material  Bodies,"  Report  TR-74-7,  Dept,  of  Electrical  and  Com¬ 
puter  Engineering,  Syracuse  University,  Syracuse,  N.Y.,  October  1974. 

[3]  T.  K.  Wu,  "Electromagnetic  Scattering  from  Arbitrarily-Shaped  Lossy 
Dielectric  Bodies,"  Ph.D.  Dissertation,  University  of  Mississippi,  1976. 

[4]  C.  Muller,  Foundations  of  the  Mathematical  Theory  of  Electromagnetic 
Waves,  Springer-Verlag,  1969,  p.  301,  Equations  (40)-(41).  (There  are 
some  sign  errors  in  these  equations.) 

[5]  V.  V.  Solodukhov  and  E.  N.  Vasil’ev,  "Diffraction  of  a  Plane  Electromagnetic 
Wave  by  a  Dielectric  Cylinder  of  Arbitrary  Cross  Section,"  Soviet  Physics  - 
Technical  Physics,  vol.  15,  No.  1,  July  1970,  pp .  32-36. 

[6]  N.  Morita,  "Analysis  of  Scattering  by  a  Dielectric  Rectangular  Cylinder  by 
Means  of  Integral  Equation  Formulation,"  Electronics  and  Communications  in 
Japan,  vol.  57-B,  No.  10,  October  1974,  pp.  72-80. 
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of  revolution  by  Vasil'ev  and  Materikova  [7].  We  will  call  this  choice 
the  Miiller  formulation.  Conditions  for  the  uniqueness  of  solutions  are 
established  in  terms  of  the  combination  constants.  It  is  found  that 
solutions  to  both  the  PMCHW  formulation  and  to  Muller's  formulation  are 
unique  at  all  frequencies. 

Numerical  solutions  to  the  coupled  pair  of  equations  are  obtained  by 
the  method  of  moments  [8].  It  is  relatively  easy  to  obtain  numerical  solu¬ 
tions  to  these  equations  because  the  required  operators  are  the  same  as  those 
evaluated  in  earlier  reports  [9,  10].  An  exemplary  computer  program  capable 
of  obtaining  the  solution  to  both  the  PMCHW  formulation  and  the  Muller  formu¬ 
lation  is  described  and  listed.  This  is  a  main  program  which  uses  subroutines 
similar  to  those  in  [10]  to  compute  the  equivalent  electric  and  magnetic 
currents  and  the  two  principal  plane  scattering  patterns  for  a  loss-free 
homogeneous  body  of  revolution  excited  by  an  axially  incident  electromagnetic 
plane  wave.  Computed  results  for  the  equivalent  currents  and  principal  plane 
scattering  patterns  of  a  dielectric  sphere  whose  relative  dielectric  constant 
is  four  show  reasonable  agreement  between  our  solution  to  the  PMCHW  formula¬ 
tion,  our  solution  to  the  Muller  formulation,  and  the  "exact"  series  [11] 
solution  in  the  resonance  region.  Computer  program  subroutines  which  calcu¬ 
late  the  "exact"  series  solution  for  perfectly  conducting  spheres  as  well  as 
for  loss-free  homogeneous  spheres  will  be  described  and  listed  in  a  subsequent 
report. 


[7]  E.  N.  Vasil'ev  and  L.  B.  Materikova,  "Excitation  of  Dielectric  Bodies  of 
Revolution,"  Soviet  Physics  -  Technical  Physics,  vol.  10,  No.  10,  April 
1966,  pp.  1401-1406. 

[8]  R.  F.  Harrington,  Field  Computation  by  Moment  Methods,  Macmillan  Co., 

New  York,  1968. 

[9]  J.  R.  Mautz  and  R.  F.  Harrington,  "H-Field,  E-Field,  and  Combined  Field 
Solutions  for  ■‘Bodies  of  Revolution,"  Interim  Technical  Report  RADC-TR- 
77-109,  Rome  Air  Development  Center,  Griffiss  Air  Force  Base,  New  York, 
March  1977. 

[10]  J.  R.  Mautz  and  R.  F.  Harringtpn,  "Computer  Programs  for  H-Field, 
E-Field,  and  Combined  Field  Solutions  for  Bodies  of  Revolution," 

Interim  Technical  Report  RADC-TR-77-215 ,  Rome  Air  Development 
Center,  Griffiss  Air  Force  Base,  New  York,  June  1977. 

[11]  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGraw-Hill 
Book,  Co.,  1961.  Section  6-9. 
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II.  SURFACE  INTEGRAL  EQUATION  FORMULATION 

An  electromagnetic  field  propagating  in  a  homogeneous  medium  of 
permeability  pg  and  permittivity  is  incident  on  the  surface  S  of  a 
homogeneous  obstacle  of  permeability  p^  and  permittivity  e^.  The  sub¬ 
script  e  denotes  exterior  medium  and  the  subscript  d  denotes  diffracting 

s  s 

medium.  We  wish  to  calculate  the  scattered  electromagnetic  field  E  ,  H 
outside  S  and  the  diffracted  electromagnetic  field  E,  H  inside  S  in  terms 
of  the  electromagnetic  field  E*,  H1  which  would  exist  on  S  in  the  absence 
of  the  obstacle.  This  original  problem  is  shown  in  Fig.  1  where  J*,  M* 
are  the  electric  and  magnetic  sources  of  E1,  H1  and  n  is  the  unit  normal 
vector  which  points  outward  from  S. 

The  equivalence  principle  (stated  in  Appendix  A)  is  used  to  piece 

s  s 

together  an  outside  situation  consisting  of  medium  p  ,  e  and  field  E  ,  H 

0  0  AV 

outside  S  and  an  inside  situation  consisting  of  medium  p  ,  e  and  field 
i  i  e  e 

-E  ,  -H  inside  S.  This  composite  situation  is  shown  in  Fig.  2.  Since 

s  s  i  i 

E  ,  H  is  source-free  outside  S  and  E  ,  H  is  source-free  inside  S,  the 

only  sources  in  Fig.  2  are  the  equivalent  electric  surface  current  J  and 

the  equivalent  magnetic  surface  current  M  on  S. 

As  a  second  application  of  the  equivalence  principle,  we  combine  an 
outside  situation  consisting  of  medium  p^,  and  zero  field  with  an  inside 
situation  consisting  of  medium  p^,  and  field  E,  EL  This  combination  of 
situations  is  shown  in  Fig.  3.  Since  E,  H  is  source-free  inside  S,  the  only 
sources  in  Fig.  3  are  the  equivalent  electric  surface  current  -\J  and  the 
equivalent  magnetic  surface  current  -M  on  S.  By  using  (A-l)  and  (A-2)  to 
express  the  surface  currents  in  terms  of  the  discontinuities  of  the  tangen¬ 
tial  fields  across  S  and  by  using 

n  x  £  =  n  x  (ES  +  E1)  (1) 

Ae  -v- 

n  x  h  =  n  x  (Hs  +  H1)  (2) 

on  S,  the  interested  reader  can  verify  that  the  surface  currents  in  Fig.  3 
are  indeed  the  negatives  of  those  in  Fig.  2.  Equations  (1)  and  (2)  are  the 
boundary  conditions  that  the  tangential  components  of  the  fields  in  the 
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original  problem  as  shown  in  Fig.  1  are  continuous  across  S. 

The  scattered  field  ES,  HS  outside  S  and  the  diffracted  field  E,  H 
inside  S  could  easily  be  calculated  if  J  and  M  were  known  because  the  media 
into  which  and  M  radiate  is  homogeneous  in  Figs.  2  and  3.  We  have  to 
determine  J  and  M.  The  equivalence  principle  states  that  there  exist  J  and 
M  which  radiate  the  fields  in  Figs.  2  and  3,  but  the  equivalence  principle 
does  not  tell  what  J  and  M  are.  The  equivalence  principle  does  state  that 


J  =  n  x  H 


(3) 


M  =  E  x  n 

A/- 

but  this  is  not  very  useful  because  E  and  H  are  unknown. 
From  Figs.  2  and  3, 

-n  x  E  =  n  x  E1 

•j*  ^ 


(4) 


(5) 


-n  x  H  =  n  x  H 
—  ~e 


-n  x  e  .  »  0 
■*-  ~d 


-n  x  h.  =  0 
o-d 


(6) 

(7) 

(8) 


where 

E  is  the  electric  field  just  inside  S  due  to  J,  M,  radiating  in  p  ,e 
-v-e  0  e 

H  is  the  magnetic  field  just  inside  S  due  to  J,  M,  radiating  in  p  ,  e 
~-e  —  ^  e  e 

E,  is  the  electric  field  just  outside  S  due  to  J,  M,  radiating  in  P,»e, 

<wQ  U.  U 

H*  is  the  magnetic  field  just  outside  S  due  to  J,  M,  radiating  in  P,*£, 

— d  ■—  d  d 

The  equivalent  currents  J,  M  which  appear  in  Figs.  2  and  3  satisfy  (5) — (8) 
because  (5) -(8)  were  obtained  from  Figs.  2  and  3.  It  is  shown  in  Appendix  B 
that  the  solution  to  (5)-(8)  is  unique.  Therefore,  (5) — ( 8)  uniquely  deter¬ 
mine  the  equivalent  currents  J,  M_ of  Figs.  2  and  3. 

Equations  (5) -(8)  form  a  set  of  four  equations  in  the  two  unknowns 
J  and  M.  The  usual  methods  of  equation  solving  apply  only  when  the  number 

A*-  AA- 

of  equations  is  equal  to  the  number  of  unknowns.  We  want  to  reduce  the 
set  of  four  equations  (5) -(8)  to  two  equations.  One  way  to  do  this  is  to 
form  the  linear  combination 
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(9) 


-n  x  (E_  +  afit)  =  n  x  E1 
•***€.  -*<1  ^ 


of  (5)  and  (7)  and  the  linear  combination 


-n  x  (H“  +  gH*)  =  n  x  H1 


of  (6)  and  (8)  where  a  and  8  are  complex  constants. 

The  solution  J,  M  to  (5) — ( 8)  satisfies  (9)  and  (10).  This  J,  M 
will  be  the  only  solution  to  the  pair  of  equations  (9)  and  (10)  if 


-n  x  (E  +  aE.) 
•**-  -w-a 


-n  x  (H  +  gH  )  =  0 

■v*-  <1 

have  only  the  trivial  solution  J  =  M  =  0.  From  (11)  and  (12), 


P  =  -  ccg  P  (13) 

e  d 

where  P  is  the  complex  power  flow  of  E  ,  H  inside  S  and  P,  is  the  complex 
power  flow  of  E^,  outside  S.  The  asterisk  in  (13)  denotes  complex  con¬ 
jugate.  If  ag*  is  real,  then  the  real  part  of  (13)  reduces  to 

Real(P  )  =  -  ag*Real(P,)  (14) 

e  d 

* 

If  ag  is  not  only  real  but  also  positive,  then 


Real(P,)  =  0 
d 

because  both  Real(Pe)  and  Real(Pd)  are  greater  than  or  equal  to  zero. 
Equation  (15)  implies  that 

n  «  E  j  =  n  x  H  j  =  0 

As- 

Substitution  of  (16)  into  (11)  and  (12)  yields 


nxE  =  n  x  H  =  0 
~e  •—  — e 
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The  system  of  equations  (16)  and  (17)  is  precisely  the  homogeneous  system 
of  equations  associated  with  (5)-(8).  It  was  shown  in  Appendix  B  that 
this  homogeneous  system  of  equations  has  only  the  trivial  solution  J  -  M_  **  0. 
Therefore,  if  a8  is  real  and  positive,  then  the  coupled  pair  of  equations 
(11)  and  (12)  has  only  the  trivial  solution  J  =  M  =  0  so  that  the  solution 
J,  M  to  (5)-(8)  is  the  only  solution  to  the  coupled  pair  of  equations  (9) 
and  (10) . 

If  a  =  8  =  1,  then  (9)  and  (10)  become 

-n  x  (E“  +  E_t)  =  n  x  e1  (18) 

-n  x  (H-  +  H^)  =  n  x  h1  (19) 

-*-0  vw- 

The  set  of  equations  (18)  and  (19)  is  the  coupled  pair  of  surface  integral 

equations  transcribed  by  Poggio  and  Miller  [1].  We  call  these  equations 

* 

the  PMCHW  equations.  Since  a  =  8  =  1  implies  that  aB  is  real  and  positive, 
the  argument  consisting  of  (11) -(17)  and  involving  real  power  flow  shows 
that  (18)  and  (19)  uniquely  determine  the  desired  J,  M  of  Figs.  2  and  3. 

That  (18)  and  (19)  uniquely  determine  J,  M  of  Figs.  2  and  3  can 
also  be  shown  as  follows.  The  desired  J,  M  of  Figs.  2  and  3  satisfies 
(18)  and  (19)  because  (18)  and  (19)  were  obtained  from  Figs.  2  and  3. 

This  desired  J,  M  will  be  the  only  solution  to  (18)  and  (19)  if  the  associ- 

\A* 

ated  set  of  homogeneous  equations 

-n  x  (e~  +  E*)  =  0  (20) 

— e  -~*a 

-n  x  (h~  +  H*)  =  0  (21) 

~e  ~*d 

has  only  the  trivial  solution  J  =  M  =  0. 

The  following  argument  shows  that  (20)  and  (21)  have  only  the  trivial 

solution  J  =  M  =  0.  Let  E,,  H,  be  the  electromagnetic  field  outside  S  due 
—  — -  «^a 

to  J,  M  radiating  in  yd>  ed-  Let  E^,  Hg  be  the  electromagnetic  field  inside 
S  due  to  J,  M  radiating  in  y  ,  eg.  Use  the  equivalence  principle  to  form 
the  composite  situation  consisting  of  medium  yd,  cd  and  field  E^,  outside 
S  and  medium  y  ,  e  and  field  -E  ,  — H  inside  S  as  shown  in  Fig.  4.  In 

6  Q.  '**-6  '“■c 

Fig.  4,  Ed>  Hd  is  a  source-free  Maxwell '.an  field  outside  S.  Since  E e,  is 
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a  source-free  Maxwellian  field  inside  S,  the  field  -E  ,  -Hg  appearing  in 
Fig.  4  is  also  a  source-free  Maxwellian  field  inside  S.  Now,  (20)  and  (21) 
state  that  the  tangential  components  of  the  field  in  Fig.  4  are  continuous 
across  S.  Thus,  Fig.  4  is  entirely  source-free  so  that  the  field  in  Fig.  4 
is  zero  everywhere  in  which  case  (B-l)-(B-4)  are  satisfied.  But,  as  shown  in 
Appendix  B,  (B-l)-(B-4)  have  only  the  trivial  solution  J  =  M_=  0.  Hence, 

(20)  and  (21)  have  only  the  trivial  solution  J  =  M  =  0. 

If 


(22) 


then  (9)  and  (10)  become 

£  j  »  J 

-n  x  (E~  -  —  E .)  =  n  x  E 
-  «~e  c  —a  ~~ 

e 


(23) 


(24) 


-n  x  (H 


—  H*)  -  n  x  H1 
p  —d  — 

e 


(25) 


The  set  of  equations  (24)  and  (25)  is  the  coupled  pair  of  surface  integral 
equations  obtained  by  Muller  [4].  We  call  these  equations  the  Muller  equa¬ 
tions.  The  singularity  that  4_he  kernels  of  the  integral  equations  (24) 
and  (25)  exhibit  as  the  source  point  passes  through  the  field  point  is  not 

as  pronounced  as  the  singularity  of  the  kernels  of  (18)  and  (19) .  If 

* 

u  ,  e  ,  v,  ,  and  c,  are  real  in  (22) -(25),  then  aB  is  real  and  positive, 
e  e  d  a 

In  this  case,  the  argument  consisting  of  (11) -(17)  shows  that  (24)  and  (25) 
uniquely  determine  the  desired  J,  M  of  Figs.  2  and  3. 

An  alternate  proof,  valid  for  lossy  media,  that  (24)  and  (25)  uniquely 
determine  the  desired  J,  M  is  presented.  This  proof  is  similar  to  the  argu- 
ment  which  used  Fig.  4  to  show  that  (18)  and  (19)  uniquely  determine  the 
desired  J,  M  and  is  as  follows.  The  desired  JT,  _M  of  Figs.  2  and  3  satisfies 
(24)  and  (25)  because  (24)  and  (25)  were  obtained  from  Figs.  2  and  3.  This 
desired  J,  M  will  be  the  only  solution  to  (24)  and  (25)  if  the  associated 
set  of  homogeneous  equations 
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(26) 


-n  x  (E~ - -  Ej)  =  0 

•—  ~e  c  — d 
e 


-n  x  (H_  -  —  Hj)  =  0 
-  ~e  p 

e 


(27) 


has  only  the  trivial  solution  J  =  M  =  0. 

■V- 

The  following  argument  shows  that  (26)  and  (27)  have  only  the  trivial 
solution  J  =  M  =  0.  Let  E^,  be  the  electromagnetic  field  outside  S  due 
to  J,  M  radiating  in  p.,  e,.  Because  the  electromagnetic  field  E,,  H,  is 

'*’•**•  a  d  «*-d  — d 

a  source-free  Maxwellian  field  in  p^,  outside  S,  the  dual  electromagnetic 

field  n  ,  H ,,  — —  E,  where 
d  -wl  n  ,  — d 

d 


■i 


(28) 


is  also  a  source-free  Maxwellian  field  outside  S. 
electromagnetic  field  inside  S  due  to  J,  M  radiating  in  p 


Let  E  ,  H  be  the 
~e  ~e 


»  E  . 

e  e 


Because 


the  electromagnetic  field  E  ,  H  is  a  source-free  Maxwellian  field  inside  S, 

^  C  ^ 

the  dual  electromagnetic  field  n  H  ,  -  E  where 

6  D  '^G 

e 


(29) 


is  also  a  source-free  Maxwellian  field  inside  S.  Use  the  equivalence 

principle  to  form  the  composite  situation  consisting  of  medium  p_,  e,  and 

1  d  d 

field  n ,  H  ,, - E ,  outside  S  and  medium  p  ,  e  and  field 

d  Ti  -  mi  e  e 

_  d 

U  e  - 

C  6  1 

-  (n  H  ,  -  E  )  inside  S  as  shown  in  Fig.  5.  Now,  (26)  and  (27) 

p,e,  e-e  n  ~e  6  ’ 

da  e 

state  that  the  tangential  components  of  the  field  in  Fig.  5  are  continuous 
across  S.  Thus,  Fig.  5  is  entirely  source-free  so  that  the  field  in  Fig.  5 
is  zero  everywhere  in  which  case  (B-l)-(B-4)  are  satisfied.  But,  as  shown 
in  Appendix  B,  (B-l)-(B-4)  have  only  the  trivial  solution  J  =  M  =  0.  Hence, 
(26)  and  (27)  have  only  the  trivial  solution  J  =  M  =  0. 
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III.  METHOD  OF  MOMENTS  SOLUTION  FOR  A  BODY  OF  REVOLUTION 


In  this  section,  a  method  of  moments  solution  to  (9)  and  (10)  is 
developed  for  a  homogeneous  loss-free  body  of  revolution.  Special  cases  of 
(9)  and  (10)  are  the  PMCHW  equations  (18)  and  (19)  and  the  Muller  equations 
(24)  and  (25). 


For  compatibility  with  equation  (40)  on  page  14  of  [9],  we  rewrite 

(9)  as 


+  ctE+) 

~d  tan 


—  E1 
n  '-tan 
e 


(30) 


where  tan  denotes  tangential  components  on  S  and  ne  is  given  by  (29).  The 
fields  on  the  left-hand  sides  of  (30)  and  (10)  are  written  as  the  sum  of 
fields  due  to  J  and  fields  due  to  M.  Advantage  is  taken  of  the  fact  that 
the  operator  which  gives  the  electric  field  due  to  a  magnetic  current  is 
the  negative  of  the  operator  which  gives  the  magnetic  field  due  to  an 
electric  current  and  that  the  operator  which  gives  the  magnetic  field  due 
to  a  magnetic  current  is  the  square  of  the  reciprocal  of  the  intrinsic 
impedance  times  the  operator  which  gives  the  electric  field  due  to  qn 
electric  current.  In  view  of  the  above  considerations,  (30)  and  (10)  be¬ 
come 


(-  —  E  (J)  +  —  H~(M)  - 
n  -e  ~  n  ~e  — 
e  e 


a 


~  E,  (J)  + 

~  r) 
e  e 


—  e* 

n  -tan 
e 


(31) 


-n  x  (H~(J)  +  ~  E  (M)  +  BHj(J)  +  -4-  E,  (M) )  =  n  x  H1  (32) 

ne  nd 

where  E_  denotes  the  operator  which  gives  the  electric  field  due  to  an 
electric  current.  The  subscript  e  or  d  on  E  denotes  radiation  in  either 
ye»  eg  or  The  superscript  +  or  -,  if  present  on  E,  denotes  field 

evaluation  either  just  outside  S  or  just  inside  S.  The  H's  in  (31)  and- 
(32)  are  the  corresponding  magnetic  field  due  to  electric  current  operators. 
We  stress  that  all  E’s  and  H^s  in  (31)  and  (32)  are,  by  definition,  operators 
which  give  electric  and  magnetic  fields  due  to  electric  currents,  even  though 
these  operators  act  on  both  electric  and  magnetic  currents  J  and  M  in  (31) 
and  (32). 
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Let 


N 


J  =  l  l  (It,  J*.  +  I9.  J9.) 
—  L  .  ,  nj  -nj  nj  ~ni 

n=_oo  j=i  j  j  J  J 


T  Q 


(33) 


oc  N 


m  =  n  I  l  (vc.  j\  +  v*  .) 

-  e  n=L_oo  nj  ~nj  nj  -nj 


(34) 


where  I1”  ,  1^  ,  Vfc , ,  and  V^,  are  coefficients  to  be  determined  and 
nj  nj  nj  n  j 


A 

~nj 

=  u^fj (t)e^n^ 

(35) 

J*. 

~nj 

=  uf,(t)ejn* 

J 

(36) 

In  (35)  and  (36),  t  is  the  arc  length  along  the  generating  curve  of  the  body 
body  of  revolution  and  <j>  is  the  longitudinal  angle.  ufc  and  u^  are  unit 
vectors  in  the  t  and  <J>  directions  respectively  such  that  u^  *  Ht  “  2.  and 
fj(t)  is  the  scalar  function  of  t  defined  on  page  10  of  [93*  The  body  of 
revolution  and  coordinate  system  are  shown  in  Fig.  6.  Substitution  of 
(33)  and  (34)  into  (31)  and  (32)  yields 


»  N 


y  y  -f  (h  (jc.)  +  ohW.))  vt.  + 

L  .  ,  )  -e  -ni  -d  ~nj  tan  nj 
n=-°°  j=l  *- 


(H"(J*  )  +  oH*(J<{\))  + 

-e  -*nj  ~d  ~nj  tan  nj 


E AJL)  E^A) 


E  (J*.)  an  E.(J*  ) 


W  A  Ei  .  V  J  .J  .  l,  ,  /  *  1  W  •  /  j 

(  _ 2.  rd-rnj— n  It  +  (_r£_riy _ £ 


nd  tan  nj 


nd  '  tan  nj 


}■-< 

j  n  -tan 


e 

(37) 


•  N  (  E  (A)  Bn  E.(A)  E  CJ^ .)  8np  E.(J*  )  a 

_  x  7  V  e  -nj  __e  ~d  *mj  .  yt  +  f-e  ~nj_  +  _®  r*  + 

n  1  •  1  (  n  n,  n,  nj  r\  n.  nH  nj 

n=-«  1=1  L  e  d  d  e  d  d 


(38) 


Define  the  inner  product  of  two  vector  functions  on  S  to  be  the 
integral  over  S  of  the  dot  product  of  these  two  vector  functions.  Because 
the  field  operators  in  (37)  and  (38)  are  the  same  as  those  considered  in  [9], 
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only  the  nth  term  of  the  sum  (37)  or  (38)  contributes  to  the  inner  product 

of  (37)  or  (38)  with  either  Jt  .  or  ^  ..  Hence,  the  inner  product  of  (37) 

•j—ni  **-ni 

with  J  i=l,2,...N,  and  J  i=l,2,...N,  successively  and  the  inner  pro- 
-ni  ‘ — nx  , 

duct  of  (38)  with  J  i=l,2,...N,  and  j'5’  i=l,2,...N,  successively  gives 

- — ni  ' — ni 

the  matrix  equation 


(Y*C  + 
ne 

«*£> 

<<t + 

(zfct  + 

ne 

and 

ne 

O 

(Zfc*  + 
ne 

and 

ne 

ztf 

nd 

(-y“- 

ne 

<d> 

(-Y1*  - 
ne 

<d> 

(Z*C  + 
ne 

and 

2„*$> 

(Z**  + 
ne 

and 

74><P 

nd 

(zH  +  ^  Z*J) 

ne  n,  nd 

d 


8n 


(ZH  +  z<hk 

ne  n,  nd 


Bn 


(Ytt:  +  BY^)  (Yt<P  +  BY^p 
ne  nd  ne  nd 


fttv  /Vt* 


M' 


(-ZtC  _  U®.  ztt}  (_zt$ - e  zt^  (Y^c  +  gY^)  (Y^^  +  BY^J) 

-J  '  ne  „  nd/  ne  nd'  ne 

d 


ne  n ,  nd 
d 


nd; 


) 

"  — 

P 

r 

4-1 

1  +<> 

n 

n 

) 

P 

P 

n 

n 

P 

p 

n 

n 

P 

p 

n 

n 

for  n  -  0,  +1,  +2, .  In  (39),  P ,  P,  P,  and  P  are  column 

—  —  n  n  n  n 

vectors  of  the  coefficients  appearing  in  (33)  and  (34).  Also, 


(39) 


(Ypq)  = 
Unf'ij 


JP  ,  •  n  x  Hr(Jq.)ds 
— ni  ~-f  ~nj 


(40) 


(zpq) 

{  nf'ij 


_1 

Pf  l 


f  JP  •  Ef(Jq.)ds 
J  <^-ni  — f  ~nj 


(41) 


yP  =  ±- 
ni  n 

e 


JP  .  •  E1  ds 
ni  — 


(42) 


ip.  = 

ni 


[  JP  .  ♦  n  x  H*  ds 
J  nx  *<*  — 


(43) 


where  p  may  be  either  t  or  <f>,  q  may  be  either  t  or  <j>,  and  f  may  be 

either  e  or  d.  If  p=q  in  (40),  it  matters  whether  the  magnetic  field 

H  (Jq  )  is  evaluated  just  outside  or  just  inside  S.  The  Y’s  without 
-f  ~nj 
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carets  in  (39)  are  given  by  the  right-hand  side  of  (40)  in  which  the  mag¬ 
netic  field  is  evaluated  just  inside  S.  The  Y's  with  carets  in  (39)  are 
given  by  the  right-hand  side  of  (40)  with  magnetic  field  evaluation  just 
outside  S. 

The  Y  and  Z  submatrices  on  the  left-hand  side  of  (39)  are  the  same 
as  in  equation  (88)  on  page  24  of  [9]  with  the  reservations  that  the  caret 
on  Y  denotes  magnetic  field  evaluation  just  outside  S,  and  the  extra  sub¬ 
script  e  or  d  denotes  radiation  in  either  p  ,  e  or  p,»  e,.  The  I  column 

e  e  d  d 

vectors  on  the  right-hand  side  of  (39)  are  the  same  as  in  equation  (88) 
on  page  24  of  [9]  whereas  the  V  column  vectors  in  (39)  are  the  same  as  the 
V’s  without  carets  in  equation  (88)  on  page  24  of  [9]. 

The  solution  V1",  V^,  I1",  and  to  the  matrix  equation  (39)  determines 
n  n  n  n 

the  equivalent  electric  and  magnetic  currents  J  and  M  according  to  (33)  and 
(34).  From  Fig.  2,  these  currents  radiate  in  p  ,  to  produce  the  scattered 
field  outside  S. 

IV.  FAR  FIELD  MEASUREMENT  AND  PLANE  WAVE  EXCITATION 

In  this  section,  measurement  vectors  are  used  to  obtain  the  far  field 
of  the  equivalent  surface  currents  J  and  M  radiating  in  p  ,  e  .  This  far 

0  g 

field  is  the  far  field  scattered  by  the  homogeneous  body  of  revolution.  For 
plane  wave  excitation,  the  composite  vector  on  the  right-hand  side  of  (39) 
is  expressed  in  terms  of  these  measurement  vectors. 

By  reciprocity, 

E®  •  II  =  f[  (J(r)  •  E(II  )  -  M(r)  *  H(II  ))ds  (44) 

—  J  J  ~  ~  —  -r  —  -  ~«rr 

S 

where  ES  is  the  far  electric  field  due  to  J  and  M,  II  is  a  receiving  elec- 
trie  dipole  at  the  far  field  measurement  point,  E(I£  )  is  the  electric  field 
due  to  II r,  and  H(IIr)  is  the  magnetic  field  due  to  11^..  Both  E(IIr)  and 
H(II  )  are  evaluated  at  point  r  on  S  where  r  is  the  point  at  which  the  dif- 
ferential  portion  of  surface  ds  is  located.  If  I.  is  tangent  to  the  radiation 
sphere. 
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(45) 


r  -jk  •  r 

E(U  )  =  - U  e  "r  ~ 

~  4irr  %-r 

r 


r  -jk  •  r 

H ( 1 8.  )  =  - (k  x  18,  )e  ~r 

~  ~r  4ivr  — r  — r 

r 


where  is  the  distance  between  the  measurement  point  and  the  origin  in 
the  vicinity  of  S.  Also,  kf  is  the  propagation  vector  of  the  plane  wave 
coming  from  18.  ,  k  is  the  propagation  constant  and  n  i3  the  intrinsic 
impedance  of  the  medium  outside  S.  To  simplify  the  notation  in  this  sec¬ 
tion,  we  have  omitted  the  subscript  e  from  all  parameters  dependent  on  the 
medium.  It  is  understood  that  all  far  field  measurement  vectors  and  plane 
wave  excitation  vectors  depend  only  on  the  external  medium  p  , 

Substitution  of  (33),  (34),  (45),  and  (46)  into  (44)  gives 


Es  = 

0  47rr 


--  I  (Rt<!>Vt  +  +  Rt0tC  +  iW), 

4-  nn  n  r>  n  n  n  n 


for  18  *  un  and 

— r  ^0 


w  intp 

l  (-RtGVl  -  +  iW)e'  r 

4  n  n  n  m  m  n  n 


IT  1"  t 

for  18.  =  u,  where  u„  and  u,  are  unit  vectors  in  the  6  and  <J>  directions 

~r  ~<|>  r  r 

respectively.  As  shown  in  Fig.  7,  0  and  <)>  are  the  angular  coordinates  of 

^  s  s 

the  receiver  location  at  which  I8r  is  placed.  In  (47)  and  (48),  E  and  E 


are  the  0  and  $  components  of  .  Also,  V* ,  V**,  Ik,  and  are  column 
r  r  —  n  n  n  n 

vectors  of  the  coefficients  appearing  in  (33)  and  (34) .  Furthermore,  R"  is 


a  row  vector  whose  jth  element  is  given  by 


~jn,t,r  f 


-Ik  ‘  r 

u  e  ds 

-q 


where  p  may  be  either  t  or  <j>  and  q  may  be  either  ft  or  (j> .  Tn  view  of  (35) 
and  (36),  (49)  is  the  same  as  equation  (92)  on  page  26  of  [9].  Tt  is 
shown  in  [9]  that  the  right-hand  side  of  (49)  does  not  depend  on  <J>  . 


«  j  "fsww  if  w  •«» 


. ,„.. ,,..  .,,,.  ^v/s^n  raw«f 


For  plane  wave  incidence  and  expansion  functions  J^ ,  and  J ^  given 

*"n  j  -nj 

by  (35)  and  (36),  the  equivalent  currents  (33)  and  (34)  and  the  fields  (47) 
and  (48)  have  special  forms.  To  obtain  these  forms,  assume  that  the  incident 
electromagnetic  field  E^,  is  either  a  0  polarized  field  defined  by 


-ik  •  r 
r,x  .  t  J*-t  -- 
L  =  kn  u  e 


ik. 


H 


=  -  k  u  e 

~y 


(50) 

(51) 


or  a  J  polarized  field  defined  by 

i 


-jk 
J  t 


li  =  kn  u  e 

~  ~y 

H  -  k  u.,  e 


(52) 


(53) 


where  k  is  the  propagation  vector  and,  as  shown  in  Fig.  7,  u.  and  u  are 
.v-t  **-y 


unit,  vectors  in  the  0  and  y  directions  respectively. 


Here,  0t  is  the 


colatitude  of  the  direction  from  which  the  incident  wave  comes. 


k^  is  in 


the  xz  plane.  No  generality  is  lost  by  putting  k  in  the  xz  plane  because 
if  kt  were  shifted  out  of  the  xz  plane  by  an  angle  <J>t»  the  response  would 
also  be  shifted  by  the  same  angle  <f>  . 

Substituting  (50)  and  (51)  into  (42)  and  (43),  then  substituting  (52) 
and  (53)  into  (42)  and  (43),  next  taking  advantage  of  the  relationships 


J  ,  x  n 
* — n  i  ~ 


=  J 


>■— n  i 


Jr  .  x  n  - 
- — nx 


Tt 

v~{  j 

^111 


(54) 

(55) 


which  are  apparent  from  (35),  (36)  and  Fig.  6,  then  comparing  the  results 
with  (49),  and  finally  using  equation  (104)  on  page  29  of  (9],  we  obtain 
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< 


“ 

*>G 

v 

'  r.B 

n 

n 

n 

n 

se 

n 

n 

n 

n 

|te 

‘ft 

-R^9 

n 

n 

n 

n 

■>t  0 
-R 

n 

n 

n 

n 

The  first  superscript  on  and  I  in  (56)  is  the.  superscript  which  appears 

on  the  right-hand  side  of  (39) .  The  second  superscript  on  and  1  in  (56) 

n  n 

denotes  the  polarization  of  the  incident  plane  wave.  If  this  second  super¬ 
script  is  6,  the  0  polarized  field  given  by  (50)  and  (51)  is  incident.  If 
this  second  subscript  is  <J>,  the  <p  polarized  field  given  by  (52)  and  (53)  is 
incident.  The  jth  element  of  the  column  vector  on  the  right-hand  side 
of  (56)  is  given  by  (49)  with  6^  replaced  by  0^,.  Conceding  that  9  does  not 
appear  explicitly  in  (49),  we  really  mean  that  0  is  replaced  by  0^  after  the 
surface  integral  in  (49)  is  evaluated.  In  other  words,  0^  is  replaced  by  0^ 
in  equation  (95)  on  page  27  of  [9]. 

For  plane  wave  incidence,  the  Fn  and  -n  terms  in  formulas  (33)  and 
(34)  for  the  equivalent  currents  can  be  combined  as  follows.  According  to 
equations  (102)  and  (103)  of  [9],  the  Y  and  Z  submatrices  in  (39)  are  either 
even  or  odd  in  n.  Hie  even-odd  properties  In  n  of  the  submatrices  of  the 
square  matrix  on  the  left-hand  side  of  (3y)  are  tabulated  as 


-  F  F  - 

F  -  —  F 

+  F 

+  -  -  F 

where  F  denotes  an  even  s  ibmatrix  and  -  denotes  an  odd  submatrix.  It 
follows  that  the  submatrices  of  the  inverse  of  the  square  matrix  on  the 
left-hand  side  of  (39)  have  even-odd  properties  In  n  given  by 
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- -•  — . 


1 


+ 


From  (56)  and  the  even-odd  properties  of  given  by  equation  (104)  or 
page  29  of  [9],  the  column  vectors  and  1  on  the  right-hand  side  of 
(39)  are  either  even  or  odd  in  n.  The  even-odd  properties  of  the  sub- 
matrices  on  the  left-hand  side  of  (56)  are  tabulated  as 


Because  of  the  above  even-odd  properties  of  tho  square  matrix  on  the  left- 
hand  side  of  (39)  and  the  column  vector  on  the  right-hand  side  of  (39), 
the  solutions  to  (39)  satisfy 


The  first  superscript  on  the  column  vectors  V  and  1,  in  (57)  is  that 

'  +T1  +TI 

which  appears  on  the  colunui  vectors  $  and  f  on  the  left-hand  side  of 
'  r  n  n 

(39).  The  second  superscript  on  the  column  vectors  in  (57)  denotes  either 
the  0  or  the  4)  polarized  in.  I  dent  plane  wave.  Substitution  of  (57),  (35), 
and  (36)  into  (33)  and  (34)  yields 


which  appears  on  the  colunui  vectors  V  and  f  on  the  left-hand  side  of 
'  r  n  n 
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>>«■—»■ .  . . . . . .  . .  .  ■ 


J8  = 

+  l 

2(flt8)u  cos(n<j>)  + 
n  —t 

n=l 

-  M9 

f)  ■-*- 

-  <f*!V  +  ? 

o  - , 

nal 

2j  (fVt8)u  sin(ntfi) 
n  -4; 

X<t>" 


n  '~4> 


(58) 


(59) 


for  the  G  polarized  incident  wave  and 


J*  "  < $!*>»*  +  f 

°  ^  n=i 

—  =  (fVt<!>)u  +  l 

n  —  '  o  -t  Ly 

n=l 


2j(f?**)u  sin(n<p)  +  2(?f**)u  cos  (nd>) 
n  v*-t  n  ~<j> 

2(fVt<^>)u  cos(n<j>)  +  2j  (fV^)u,sin(n$) 
n  ~a.  ~  n  ~4> 


(60) 


(61) 


for  the  <j>  polarized  incident  wave.  In  (58)-(61),  f  is  a  row  vector  of 
the  f  (t)  . '  The  superscript  8  or  cfi  on  J  or  M  in  (58)-(61)  differentiates 
the  equivalent  currents  for  the  9  polarized  incident  wave  from  those  for 
the  <p  polarized  incident  wave. 

The  far  scattered  fields  (47)  and  (48)  are  specialized  to  the  6 
polarized  incident  plane  wave  by  appending  the  additional  subscript  0  to 
ES  on  the  left-hand  sides  of  (47)  and  (48)  and  the  additional  superscript 
G  to  V^,  ft  and  on  the  right-hand  sides  of  (47)  and  (48).  More- 
over,  in  view  of  equation  (104)  on  page  29  of  [9]  and  (57),  the  +n  and  -n 
terms  in  (47)  and  (48)  can  be  combined.  As  a  result,  (47)  and  (48)  become 

-jkr 

E8  »  rJiF: _ l.\  +  RtGf!()+  2  y  (Rt<J^tO+  R<K$<J>0  + 

00  4ii  r  loo  o  o  *■  n  n  n  n 

r  n-1 


+ 


RC"ftM4  it*0*'*6 
n  n  n  n 


)  cos  (rt<f» 


r)| 


(62) 


i  k  r 


Jle... _ 

2  ?f  r 


y 

n=l 


-iW' 

n  n 


f  K 


^ftO 
n 
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for  the  0  polarized  Incident  plane  wave.  Similarly,  (47)  and  (48) 
become 
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for  the  <p  polarized  incident  plane  wave.  The  first  subscript  on  E  on 
the  left-hand  sides  of  (62) -(65)  denotes  the  receiver  polarization  and 
the  second  subscript  on  E  denotes  the  transmitter  polarization. 

The  scattering  cross  section  o  is  defined  by 

pq 

/  2 1 _s  |  2 

4rr  E 
r1  pq' 


g 

where  p  is  either  0  or  <p  and  q  is  either  0  or  $.  In  (66),  E  is  e  com- 

i  ii  ^ 

ponent  of  the  scattered  field  given  by  (62) -(65)  and  j E  |  is  the  magnitude 
of  the  electric  field  of  the  incident  plane  wave.  According  to  (50)  and  (52), 


E1  =  kri 


for  both  polarizations  so  that 


4rrr2|ES  [2 
r1  pq' 


Normalized  versions  of  (68)  are 


a  4r2  E  2 

_es.  =  r  pq  . 
2  ,222 
ira  k  a  n 


a  r2 j ES  |2 

m  -*1* S-  (70) 

X  TTT) 

where  a  is  some  characteristic  length  associated  with  the  scatterer. and  X  is 
the  wavelength  in  the  external  medium. 


V.  EXAMPLES 


A  computer  program  has  been  written  to  calculate  the  equivalent 
currents  and  scattering  patterns  for  a  dielectric  body  of  revolution  ex¬ 
cited  by  an  axially  incident  plane  wave.  This  program  is  described  and 
listed  in  Part  Two.  Some  computational  results  obtained  with  this  program 
are  given  in  this  section. 


Figures  8  and  9  show  the  magnitude  and  phase  of  the  normalized 

Je  J<t  Me  Mi> 

equivalent  currents  —  ,  rf-  ,  — ,  and  —  on  the  surface  of  a  dielectric 

H  H  &  h 

y  y  x  x 

sphere  for  which  ka  =  3  and  =  4.  Here,  k  is  the  propagation  constant 

in  free  space,  a  is  the  radius  of  the  sphere  and  e  is  the  relative 

dielectric  constant  of  the  sphere.  Figure  8  represents  our  solution  of  the 
PMCHW  formulation.  Figure  9  depicts  our  solution  of  the  Muller  formulation. 


In  Figs.  8  and  9,  the  incident  field  is  a  plane  wave  traveling  in 
the  positive  z  direction.  THETA  =  0°  is  the  forward  scattering  direction 
and  THETA  =  180°  is  the  backs cat tering  direction.  The  incident  field  is 
given  by  (50)  and  (51)  with  0  =  180°.  The  origin  r  =  0  is  at  the  center 

of  the  sphere.  In  Figs.  8  and  9,  Jg  is  the  uQ  =  -  ut  component  of  electric 

current  (58)  versus  0  in  the  <fc  =  0  plane,  J,  is  the  u,  component  of  (58) 

m 

versus  0  in  the  <p  =  90°  plane,  M  is  the  u0  =  -  u  component  of  magnetic 

current  (59)  versus  0  in  the  <f>  =  90°  plane,  and  M^  is  the  u^  component  of 

(59)  versus  0  in  the  4>  =  0°  plane.  For  axial  incidence,  only  the  n=l  term 
is  present  in  (58)  and  (59).  The  symbols  x  and  +  denote  respectively 
magnitude  and  phase  of  the  method  of  moments  solution  for  the  pertinent 
component  of  the  electric  or  magnetic  current.  The  solid  curves  are  the 
exact  equivalent  currents  obtained  from  the  Hie  series  solution  [11].  The 
normalizing  constants  E^  and  H  are  defined  in  terms  of  the  incident  field 
(50)  and  (51)  by 


E 

x 


r=0 


H  =  u 

y 


r=0 


k 


(71) 


[11]  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGraw-Hill 
Book  Co.,  1961.  Section  6-9. 
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where  k  and  n  are  respectively  the  propagation  constant  of  free  space  and 
the  impedance  of  free  space. 

The  currents  of  Figs.  8  and  9  were  obtained  by  using  a  20  point 
Gaussian  quadrature  formula  for  all  integrations  in  4>.  All  Integrations 
over  the  functions  {f  (t)}  in  t  =  a(-n-9)  were  done  by  sampling  each 
fj(t)  four  times.  The  (f^(t)}  consisted  of  14  overlapping  triangle 
(divided  by  the  cylindrical  coordinate  radius)  functions  equally  spaced 
in  8.  More  precisely, 

NP  =  31 

NPHI  «  20  (72) 

MT  =  2 

where  the  above  variables  are  input  data  for  the  computer  program  described 
and  listed  in  Part  Two,  Section  V. 

Figures  10  and  11  show  the  scattering  patterns  radiated  by  the  cur¬ 
rents  of  Figs.  8  and  9  respectively.  The  symbols  X  and  +  denote 

O0Q  ^(f)  0 

— j  anc*  "2  respectively.  The  solid  curves  are  the  exact  patterns  obtained 
xa  Tra 

from  the  Mie  series  solution  [11  jj.  The  patterns  a.,  and  a,n  are  given  by 

u  U  <po 

(69),  (62)  and  (63).  Here  a  is  the  0  polarized  pattern  versus  0  in  the 

u  v  T 

<}j  =  0  plane  and  a.,  is  the  <)>  polarized  pattern  versus  0  in  the  "  90° 

0(7  r 

plane.  Tlie  THETA  In  Figs.  10  and  11  refers  to  0^.  For  axial  incidence, 
only  the  n=l  terms  are  present  in  (62)  and  (63).  Elsewhere  [3,  12],  the 
pattern  is  called  the  horizontal  polarization  because  it  is  polarized 
parallel  to  the  scattering  plane.  Similarly,  the  pattern  o'  is  called  the 
vertical  polarization  because  it  is  polarized  perpendicular  to  the  scatter¬ 
ing  plane. 

Figures  12-17  show  the  scattering  patterns  for  three  other  dielectric 
spheres.  Figures  12  and  13  are  for  relative  dielectric  constant  e  “  1.1, 
Figs.  14  and  15  for  =  10.,  and  Figs.  16  and  17  for  -  20.  All  other 

[12]  P.  Barber  and  C.  Yeh,  "Scattering  of  Electromagnetic  Waves  by 

Arbitrarily  Shaped  Dielectric  Bodies,"  Applied  Optics,  vol.  14,  No.  12, 
December  1975,  pp.  2864-2872, 
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Fig.  10.  Plane  wave  scattf  ring  patterns 
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Fig,  .12.  Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  =  3,  t  «  1.1, 
PMCHW  soJutlcn'i.  Symbols  *  and  +  denote  horizontal  polarization  and 
vertical  polarization  respectively.  Solid  line  denotes  exact  solution. 
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Fig,  13.  Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  *  3,  e  *  1.1 
Muller  solution.  Symbols  y  and  +  denote  horizontal  polarization  and 
vertical  polarization  respectively.  Solid  line  denotes  exact  solution 
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tering  patterns  for  dielectric  sphere,  ka  »  3,  ,  -  20, 

Symbols  '  and  +  denote  horizontal  pol arizat ion rand 
zat ion  respectively.  Solid  line  denotes  exact  solution 
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Plane  wave  smattering  patterns  for  dielectric  sphere,  ka  -  3,  r ^  “  20, 
Mil  Her  solution.  Symbols  •  and  +  denote  horizontal  polarization  and 
vertical  polarization  respectively.  Solid  line  denotes  exact  solution. 
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parameters  in  Figs.  12-17  are  the  same  as  in  Figs.  10  and  11.  In  Figs.  12 
and  13,  values  less  than  0.0001  are  plotted  at  0.0001. 

Figures  18  and  19  show  the  computed  scattering  patterns  of  a  finite 
dielectric  cylinder  of  radius  a  and  height  2a  when  a  is  0.25  free  space 
wavelengths.  The  relative  dielectric  constant  of  the  cylinder  is  =  4. 

The  incident  field  is  a  plane  wave  traveling  in  the  positive  z  direction, 
the  same  field  which  was  incident  upon  the  previous  dielectric  spheres. 

^00  OA0 

Figure  18  shows  — ^  an^  ~  2  as  °^ta^ne<^  from  our  solution  of  the  PMCHW 
ita  ira  g 

formulation.  Figure  19  shows  — and  as  obtained  from  our  solution 

tt a  1:3  a  a 

of  the  Muller  formulation.  The  patterns  — —■  and  are  plotted  with  the 

■rra  ira 

symbols  *  and  4-  respectively. 

The  equivalent  currents  which  radiate  the  patterns  of  Figs.  18  and  19 
were  obtained  by  using  a  48  point  Gaussian  quadrature  formula  for  all  inte¬ 
grations  in  <}>.  All  integrations  in  t  over  the  functions  (f^(t)}  were  done  by 
sampling  each  f ^ (t)  four  times.  The  { f ^  C t) >  consisted  of  11  overlapping 
triangle  (divided  by  the  cylindrical  coordinate  radius)  functions  equally 
spaced  in  t.  More  precisely, 

NP  =  25 

NPHI  =  48  (73) 

MT  =  2 

where  the  above  variables  are  input  data  for  the  computer  program  described 
and  listed  in  Part  Two,  Section  V. 

VI.  DISCUSSION 

According  to  Figs.  12  and  13,  -he  scattering  patterns  obtained  from 
our  solution  of  the  MUller  formulation  are  more  accurate  than  those  ob~  . 
tained  from  our  solution  of  the  PMCHW  formulation  for  the  dielectric  sphere 
with  ka  =  3  and  e  =1.1.  From  plots  not  included  in  this  report,  we  ob¬ 
served  that  both  our  PMCHW  solution  and  our  Muller  solution  for  the  equivalent 
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currents  on  the  dielectric  sphere  were  reasonably  accurate.  However,  the 

following  argument  shows  that  when  is  near  one,  a  slight  inaccuracy  in 

the  equivalent  currents  could  affect  the  scattering  patterns  drastically. 

As  e  approaches  one,  the  equivalent  electric  and  magnetic  currents 
r  i  i 

approach  n  *  H  and  E  x  n  respectively  whereas  the  scatterning  patterns 
approach  zero.  This  means  that  the  equivalent  electric  and  magnetic 
currents  produce  fields  which  nearly  cancel  each  other.  Hence,  a  slight 
inaccuracy  in  the  equivalent  currents  could  cause  a  large  percentage 
inaccuracy  in  the  scattering  patterns. 

We  believe  that  our  Muller  solution  is  more  accurate  than  our 
PMCHW  solution  whenever  is  close  to  one.  When  a  and  0  are  given  by 
(22)  and  (23)  as  in  the  Muller  formulation,  the  left-hand  sides  of  (9) 
and  (10)  approach  -M  and  respectively  as  approaches  one.  In  this 
case,  the  expected  solution 

J  =  n  *  H1 

M  =  E"-  x  n 

^  AKr 

can  be  obtained  by  inspection  of  (9)  and  (10).  However,  if  a  -  B  =  1 
as  in  the  PMCHW  formulation,  the  solution  to  (9)  and  (10)  is  not  obvious 
when  .  e  =  1  because  the  field  operators  on  the  left-hand  sides  of  (9) 
and  (10)  are  not  diagonal.  With  our  Muller  solution,  the  matrix  on  the 
left-hand  side  of  (39)  would  become  tridiagonal  for  t  =  1  if  its  first 
two  rows  of  submatrices  were  interchanged.  With  our  PMCHW  solution,  no 
such  simplification  of  this  matrix  is  possible  for  =  1. 

We  recommend  at  least  10  expansion  functions  per  wavelength  per 
component  along  the  generating  curve  of  the  dielectric  body  of  revolution. 
For  example,  if  the  generating  curve  were  one  wavelength  long,  the  order 
of  the  square  matrix  on  the  left-hand  side  of  (39)  should  be  at  least  36. 
The  number  36  is  arrived  at  as  follows.  There  should  be  at  least  9  expan¬ 
sion  functions  per  component  of  current.  We  say  9  expansion  functions 
rather  than  10  because  we  are  using  overlapping  triangle  functions  with  no 
peak  of  triangle  function  at  either  ends  of  the  generating  curve.  There 
are  two  components  of  electric  current  and  two  components  of  magnetic 
current. 
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According  to  equations  (20)-(23)  of  [9]  and  (58)~(61)  of  [9],  each 

element  of  the  square  matrix  in  (39)  is  a  triple  integral  consisting  of 

one  integration  with  respect  to  <f>  and  two  integrations  with  respect  to  t. 

The  <j>  integral  is  evaluated  by  using  a  Gaussian  quadrature  formula.  Each 

t  integration  is  done  by  crude  sampling  akin  to  the  trapezoid  rule.  In 

any  case,  there  should  be  at  least  10  sample  points  per  wavelength  in  the 

media  in  question.  For  instance,  if  p  is  the  largest  cylindrical 

coordinate  radius  of  the  dielectric  body  of  revolution  and  t  is  the 

max 

length  of  the  generating  curve,  then 


lOt 


NP  > 


max 


X 

lOTTp 


+  1 


(74) 


NPHI  > 


max 


-  X 


where  NP  and  NPHI  are  input  arguments  of  the  subroutine  YZ  described  and 
listed  in  Part  Two,  Section  II  and  X  is  the  wavelength  in  the  media  in 
question.  If  f  =  e  in  (40) -(41)  then  X  is  the  wavelength  in  the  external 
media,  but  if  f  -  d  in  (40) -(41)  then  X  is  the  wavelength  inside  the  dif¬ 
fracting  body  of  revolution.  The  main  program  in  Part  Two,  Section  V  is 
oversimplified  in  that  it  uses  the  same  values  of  NP  and  NPHI  for  both 
f  =  e  and  f  =  d. 

Loss  of  accuracy  in  the  computed  patterns  of  Figs.  16  and  17  may 
be  due  to  the  fact  that  (74)  was  violated.  According  to  (74),  the  values 
of  NP  and  NPHI  for  f  =  d  should  be  nearly  70  or  greater  instead  of  the  low 
values  appearing  in  (72)  .  Unfortunately,  increasing  the  values  of  the 
variables  NP  and  NPHI  increases  the  computer  time  required  to  solve  the 
problem. 

We  have  been  trying  to  obtain  accurate  numerical  results  for  the 
dielectric  sphere  for  which  a  =  0.2  free  space  wavelengths  and  =  80. 
from  our  general  dielectric  body  of  revolution  program.  We  have  not  been 
able  to  obtain  clear-cut  convergence  with  respect  to  the  variables  on  the 
left-hand  sides  of  (72)  because  we  could  not  afford  to  increase  them  as 
much  as  desired.  Our  PMCHW  solution  and  our  Muller  solution  for  the 
equivalent  currents  and  scattering  patterns  differ  from  each  other  and 
from  the  exact  solution. 
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For  the  sphere,  each  element  of  the  square  matrix  or  the  left-hand 
side  of  (39)  can  be  written  as  a  sum  over  the  infinite  set  of  spherical 
modes.  So  far,  we  have  not  been  able  to  successfully  implement  this  alter¬ 
nate  evaluation  of  the  matrix  elements  in  terms  of  spherical  modes.  The 
major  diffi<  ulty  seems  to  be  lack  of  agreement  of  a  few  matrix  elements  for 
which  both  expansion  and  testing  functions  are  near  one  of  the  poles  of  the 
sphere. 

Both  the  PMCIIW  solution  and  the  Muller  solution  are  obtained  by  taking 
a  linear  combination  of  (5)  and  (7)  and  a  linear  combination  of  (6)  and  (8) . 
There  are  two  other  possibilities  which  are 

(1)  A  linear  combination  of  (5)  and  (6)  and  a  linear  combi¬ 
nation  of  (7)  and  (8) . 

(2)  A  linear  combination  of  (5)  and  (8)  and  a  linear  combi¬ 
nation  of.  (6)  and  (/). 

These  other  two  possibilities  give  rise  to  alternative  numerical  solutions 
which  may  compare  favorably  with  the  PMC.HW  solution  and  the  Muller  solution. 
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APPENDIX  A 


THE  EQUIVALENCE  PRINCIPLE 


Let  B  ,  H  he  an  electromagnetic  field  defined  outside  a  closed 
surface  S.  The  permeability ,  permittivity,  and  electric  and  magnetic 
source  currents  outside  S  are  p  ,  i:  ,  J  ,  and  M  respectively.  This 
outside  situation  where  the  subscript  e  stands  for  "exterior  medium" 
is  shown  in  Fig.  A-l.  In  Fig.  A-l,  the  media  and  sources  inside  S  are 
undisclosed.  Let  E  ,  H  be  an  electromagnetic  field  defined  inside  S 
where  the  permeability,  permittivity,  and  electric  and  magnetic  sources 
are  p  ,,  e.,  J,,  and  M.  respectively.  This  inside  situation  where  the 
subscript  d  stands  for  "diffracting  medium"  is  shown  in  Fig.  A-2.  In 
Fig,  A-2,  the  media  and  sources  outside  S  are  undisclosed.  The  equi¬ 
valence  principle  states  that  the  solution  to  the  composite  radiation 

problem  consisting  of  medium  p  ,  t  and  sources  J  ,  M  outside  S,  medium 

e  e 

*Jd’  ed  an(*  sources  insi^e  S,  and  electric  and  magnetic  surface 

currents  J,  M  on  S  given  by 


J  =  n  x  (i!  -  H.)  (A-l) 

^-<1 

M  -  (E  -  R)  x  n  (A-2) 

where  n  is  the  exterior  unit  normal  vector  on  S  is  the  composite  elec- 

<4*. 

tromagnetie  field  E,  H  defined  by 


E  ,  H 
~e  +-e 

outside  S 

(A-3) 

V  2d 

inside  S 

(A-4) 

The  composite  radiation  problem  is  shown  in  Fig.  A- 3  which  is  entitled 
composite  situation. 

The  equivalence  principle  is  proved  by  showing  that  the  configu¬ 
ration  of  media  and  sources  in  Fig.  A- 3  gives  rise  t:o  the  composite  field 
E,  u  defined  by  (A-3)  and  (A-4) .  Now,  E,  H_  wi  1 1  be  the  field  generated 
by  the  media  and  sources  of  Fig.  A-3  if  E,  H  satisfies  Maxwell's  equations 
with  source  terms  included  and  the' radiation  condition  at  infinity.  E,  H 
satisfies  Maxwell's  equations  outside  S  and  inside  S  because  E  ,  H  and 
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Fig.  A-l.  Outside  Situation. 
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Fig.  A-2.  Inside  Situation. 
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E^,  Hj  satisfy  Maxwell’s  equations  outside  S  and  inside  S  respectively. 

E,  H  also  satisfies  the  radiation  condition  at  infinity  because  E  ,  H 
•-»  J  ~>e 

satisfies  the  radiation  condition  at  infinity.  It  remains  to  be  shown 
that  Maxwell's  equations  for  E,  exhibit  the  surface  current  sources 
J  and  M  given  by  (A-l)  and  (A-2) . 


It  is  well  known  that  a  surface  current  source  on  S  gives  rise  to 
a  discontinuity  across  S  ol:  the  tangential  component  of  the  field.  The 
preceding  statement  is  easily  verified  by  means  of  an  argument  based  on 
the  integral  forms  of  Maxwell's  equations.  Now,  this  same  argument  can 
be  construed  to  imply  that  a  discontinuity  across  S  of  the  tangential 
component  of  the  field  gives  rise  to  a  surface  current  source  on  S. 

Hence  Maxwell's  equations  for  E,  H  exhibit  the  electric  and  magnetic 
surface  current  sources  J.  and  H  on  S  given  by  (A-l)  and  (A-2).  There¬ 
fore,  E,  H  is  the  solution  to  the  composite  radiation  problem  shown  in 
Fig.  A- 3  because  E,  H  satisfies  Maxwell's  equations  with  sources  and  the 

°  /VK 

radiation  condition  at  infinity. 


APPENDIX  B 


PROOF  THAT  THE  SOLUTION  TO  EQUATIONS  (5)- (8)  IS  UNIQUE 


The  solution  J,  M  to  (!i)-(8)  will  be  unique  if  the  associated  set 
of  homogeneous  equations 


-n  x 

E  =0 
-e 

(B-l) 

-n  x 

H"  -  0 

(B-2) 

-n  x 

£-° 

(B-3) 

-n  x 

A— 

4  - 0 

(B-4) 

has  only  the  trivial  solution  J  =  M  =  0. 

vv  ^ 

From  (B-l)  and  (B~2),  J,  M  radiate  in  y  ,  e  to  produce  a  field 

"  •v'  e  e 

whose  tangential  components  are  zero  just  inside  S,  Hence,  according 

to  the  relation  between  J,  M  and  the  discontinuity  of  tangential  field 

across  S  as  exemplified  by  (A~l)  and  (A-2),  the  field  E  ,  H  radiated  by 

J,  M  in  u  ,  e  outside  S  satisfies 
~  ~  e  e 


n  x  H  -  J  (B-5) 

x  n  =  M  (B-6) 


just  outside  S.  See  Fig,  B-l. 

From  (B-3)  and  (B~4)  the  electric  and  magnetic  currents  -J,  -M  radiate 

in  y^,  r.j  to  produce  a  field  whose  tangential  components  are  zero  just  outside 

S.  Hence,  according  to  the  relation  between  -J,  -M  and  the  discontinuity  of 

tangential  field  across  S,  the  field  -E.,  -H ,  radiated  by  -J,  -M  in  y.,  e, 

-■d  d  d 

satisfies 

n  x  (-H  )  -  j  (B-7) 


( — E  , )  x  n  -  M  (B-8) 

just,  inside  S.  See  Fig.  B-2 . 

The  equivalence  principle  is  used  to  combine  the  outside  situation 
in  Fig.  B--1  with  the  inside  situation  in  Fig.  B-2  to  obtain  the  composite 
situation  shown  in  Fig.  B-3.  Because  of  (B-5)-(B-8),  the  composite 


Fig.  B-l.  Radiation  of  J,M  according  to  (B-l)  and  (B-2) . 


pig.  B-2.  Radiation  of  -J,  -M  according  to  (B-3)  and  (B-4) 


Fig.  B-3.  Composite  Situation. 


situation  In  Fig.  B-3  is  source-free •  Therefore,  the  field  in  Fig.  B-3 
is  zero  everywhere.  Hence,  the  fields  in  Figs,  B-l  and  B-2  are  zero 
everywhere  which  implies  that  J  =  =  0. 

Thus,  the  solution  to  (5) -(8)  is  unique  because  the  associated  set 
(B-l)-(B--4)  of  homogeneous  equations  has  only  the  trivial  solution. 
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APPENDIX  C 


MORE  EXAMPLES 


The  equivalent  currents  and  scattering  patterns  for  the  dielectric 
spheres  for  which  =  4  and  ka  =  4,  5,  and  6  are  plotted  In  Appendix  C. 
Figures  C-l  to  C-4  are  for  ka  =4,  Figs.  C~5  to  C-8  are  for  ka  «  5,  and 
Figs.  C-9  to  C-1.2  are  for  ka  =  6.  All  other  parameters  are  the  same  as  in 
Figs.  8-11  in  Section  V.  In  particular,  the  input  data  for  the  computer 
program  which  generated  the  method  of  moments  results  plotted  in  Figs.  C-l 
to  C-12  is  given  by  (72). 

It  is  evident  from  Figs.  8-11  and  Figs.  C-l  to  C-12  that  the  method 
of  moments  solutions  for  the  equivalent  currents  and  scattering  patterns 
are  not  as  accurate  at  ka  =  4,  5,  and  6  as  at  ka  =  3.  Loss  of  accuracy 
at  the  higher  values  of  ka  may  be  due  to  violation  of  (73). 
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Fig.  C-4.  Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  =  4,  =  4, 

Muller  solution.  Symbols  *  and  4-  denote  horizontal  polarization  and 
vertical  polarization  respectively.  Solid  line  denotes  exact  solution. 
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Fig*  C-5.  Equivalent  electric  and  magnetic  currents  for  dielectric  sphere,  ka  =5,  er  =  4,  PMCHW  solution. 
Symbols  x  and  +  denote  magnitude  and  phase  respectively.  Solid  line  denotes  exact  solution. 
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Plane  wave  scattering  patterns  for  dielectric  sphere,  ka  =•  5,  r  •  4, 
Mailer  solution.  Symbols  *  ami  +  denote  horizontal  polarization  .and 
vertical  polarization  respectively.  Sclid  line  denotes  exact  solution 
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Equivalent  electric  and  magnetic  currents  for  dielectric  sphere,  ka  =  6,  er  =  4,  M511er  solution 
Symbols  x  and  +  denote  magnitude  and  phase  respectively.  Solid  line  denotes  exact  solution. 
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PART  TWO 


COMPUTER  PROGRAM 

I.  INTRODUCTION 

The  computer  program  calculates  the  equivalent  electric  and  magnetic 
currents  (58)  and  (59)  and  the  scattering  patterns  (70)  for  a  loss-free 
homogeneous  dielectric  and/or  magnetic  material  body  of  revolution  immersed 
in  an  axially  incident  plane  wave.  This  computer  program  consists  of  a 
main  program  and  the  subroutines  YZ,  PLANE,  DECOMP  and  SOLVE. 

Part  Two  consists  of  definitions  of  the  input  and  output  for  the 
subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE,  listings  of  these  subroutines, 
definitions  of  the  input  and  output  for  the  main  program,  a  verbal  flow 
chart  of  the  main  program,  and  a  listing  of  the  main  program  with  sample 
input  and  output.  The  subroutines  YZ  and  PLANE  are  similar  to  subroutines 
of  the  same  name  in  [10].  The  subroutines  DECOMP  and  SOLVE  are,  except 
for  dimension  statements,  exactly  the  same  as  in  [13].  Hence,  the  insides 
of  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE  are  not  described  in  detail 
in  Part  Two.  Because  these  subroutines  are  quite  complicated,  a  black  box 
approach  is  suggested  wherein  the  user  is  concerned  with  just  the  input  and 
output  of  these  subroutines.  However,  the  user  is  encouraged  to  delve 
inside  the  main  program  and  to  make  any  changes  therein  that  he  deems  neces¬ 
sary  to  suit  his  needs. 


II.  THE  SUBROUTINE  YZ 
Description: 

The  subroutine  YZ (NN,  NP,  NPHI,  M,  MT,  RH,  ZH,  X,  A,  Y,  Z)  stores  the 

matrices  Y  and  Z  defined  by 
nf  nf 


nf 


Jnf 


£ 

z“ 

nf 

zt't> 

nf 

z*t 

nf 

nf 

(75) 


(76' 
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by  columns  in  Y  and  Z  respectively.  The  submatrices  on  the  right-hand 
sides  of  (75)  and  (76)  are  given  by  (40)  and  (41).  The  first  9  argu¬ 
ments  of  YZ  are  input  variables.  Except  for  the  new  input  variaoles  M 
and  MT,  the  subroutine  YZ  is  the  same  as  the  old  subroutine  YZ  on  pages 
17-21  of  [10].  If  M  *  -  1  and  MT  =  2,  these  subroutines  are  exactly  the 
same  as  far  as  the  calculation  of  Y  and  Z  in  terms  of  the  rest  of  the 
input  variables  io  concerned. 


M  ■  -  1  for  field  evaluation  just  inside  S  and  M  =  +  1  for  field 
evaluation  just  outside  5.  M  =  -  1  if  f  =  e  in  (75)  and  M  *  +  1  if 
f  “  d  in  (75).  The  value  of  M  is  not  used  in  calculating  (76)  because  the 
tangential  components  of  the  electric  field  operator  E  in  (41)  are  con-- 

,VT 


tinuous  across  S. 


All  numerical  integrations  over  t  of  f  (t)  appearing 


J.p  (35)  and  (36)  are  done  by  sampling  each  f.(t)  2*viT  times.  The  repre- 

d 

sensations  of  pf^(t)  and  —  (pf^(t))  given  b>  (6b)  Ctrl  (67)  of  [9]  ere 
replaced  by  representations  which  contain  2*MT  impulse  functions  ’nstead  of 
4  impulse  functions.  For  instance,  (66)  of  [9]  is  replaced  by 


2*MT 


pfi(t)  k  Tp+(i-l)  *2*MT  6<'t  tp-f(i-l)^MT) » 


(77) 


i-l,2,...N  where  N  will  be  defined  in  the  paragraph  which  follows  the  next 
paragraph.  The  T's  appearing  in  (77)  will  be  defined  by  (78). 

Seven  of  the  input  variables  are  the  same  an  in  the  old  subroutine 
YZ  on  pages  17-2.1  of  [10].  These  variables  are  defined  in  terms  of  variables 
appearing  in  [9]  by 


NN  ■ 

-  11 

NP  - 

P  a  ge 

9 

nth  i  ^ 

=  V 

r 

page 

1  J 

RH(l) 

;  kp  j  , 

page 

9 

ZH ( i )  = 

'  kV 

page 

9 

X(k) 

V 

page 

1  3 

A  (k)  = 

V 

page 

1  S 

60 


In  summary,  n  denotes  the  dependence  in  (35)  and  (36),  (p”,  z~) , 

i*l,2,...P  are  coordinates  on  the  generating  curve,  the  k  which  multi¬ 
plies  and  la  the  propagation  constant,  and  x^  and  are  the 
abscissas  and  weights  for  the  point  Gaussian  quadrature  integration 
in  Note  that  NP-1  should  be  an  integer  multiple  of  MT.  If  NP-1  is 

not  an  integer  multiple  of  MT,  the  program  will  ignore  RH(m)  and  ZH(m) 
for 

[(NP  -1)/MT]*MT+1  <  m  £  NP 

where  [(NP-1) /MT]  is  the  largest  integer  which  does  not  exceed  (NP-1) /MT . 

Minimum  allocations  are  given  by 
COMPLEX  Y (4*N*N) ,  Z(4*N*N) 

DIMENSION  RH(NP) ,  ZH(MP),  X(NPHI),  A(NPH1),  D(NG) 

PD(NG),  TP(2*MT*N),  CR(NPHI),  Cl(NPHI),  C2(NPHI), 

C3 (NFHI) ,  C4(NPHI) 

COMMON  RS(NG),  ZS(NG),  SV(NG) ,  CV(NG),  T(2*MT*N) 

where 

N  =•  [  (NP-1)  /MT]  -  1 
NG  =*  (N+l)  *MT 

The  variable-!  in  common  make  the  results  of  some  intermediate  calculations 
done  in  YZ  available  to  the  subroutine  PLANE  described  in  Section  III  of 
Part  Two. 

We  mention  a  few  portions  of  YZ  which  differ  from  the  subroutine 
listed  on  pages  18-21  of  [10].  Equation  (291  of  [9]  lias  been  generalized  to 


X 


2*MT*(J-1)+1  "  (  ^  dMT*(J  -l)+q  ~  2  dMT*(J-l)+I)dMT*(.I-l)+'I  '' 

k 


2  I 


1 


T 


I 


‘2*MT*J-MT-H  A.,  *A2  k  dMT*.!+q  +  2  (1MT*r+f '  dMT*  1+ 

A  q~  i. 

for  J  “  1,2,...N  and  1  "  1,2,..  ..MT  where 

MT 

1  ^  .i'i  *  (,J- 1 )  + 1 

Sfj* 


(79) 


(80) 


A2  “  k  l  dMTMti 

k  “*  1 
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(81) 


and  the  d's  are  given  by  equation  (28)  of  [9].  Note  that  A^  is  the 
electrical  length  of  generating  curve  over  which  the  first  half  of 
fj(t)  exists  and  that  A^  is  the  electrical  length  of  generating  curve 
over  which  the  second  half  of  f^.(t)  exists.  The  generalization  of 
equation  (68)  of  [9]  is 


T ' 

2*MT*(J-1)+I 


T* 

2*MT*J-MT+I 


kdMT*(J-l)+I 


'kdHT*J+I 


for  J  =  1,2,...N,  and  I  =  1,2,...MT  where  and  A^  are  given  by  (80) 
and  (81).  Expressions  (78)  -  (83)  are  calculated  in  DO  loop  68.  DO 
loop  12  accumulates  A  in  DEL.  DO  loop  19  puts  (78)  in  T(2*MT*(J-1)+I) 
and  (82)  in  TP(2*MT*(J-1)+I) .  DO  loop  15  accumulates  A2  in  DEL.  DO 
loop  16  puts  (79)  in  T(?.*MT*J-MT+I)  and  (83)  in  TP  (2*MT*J-HT+I) . 

The  subscripts  KT,  LT,  and  J1  inside  DO  loop  32  are  obtained  as 

follows.  Since  the  generating  curve  consists  of  NG  =  (N+1)*MT  small 

intervals,  it  is  composed  of  (N+l)  large  intervals  where  the  mth  large 

interval  consists  of  the  ( (m-1) *MT+1) th  through  the  (m*MT)th  small 

intervals.  The  index  I  of  DO  loop  60  denotes  the  Ith  small  interval. 

The  Ith  small  interval  is  contained  in  the  (19+1) th  large  interval  where 

19  »  ((l-D/MT].  The  second  naif  of  fjg(t)  and  the  first  half  of 

are  in  this  large  interval.  The  index  K  of  DO  loop  32  denotes  f^+K_j(t). 

Since  T( (m-1) *2*MT+1)  through  T(m*2*MT)  is  allotted  to  f  (t) ,  ra=l,2, . . . *£, 

f  (t)  is  preceded  by  (m-1)  overlaps.  For  each  overlap  the  subscript  of  T 
m 

increases  by  an  amount  MT  not  accounted  for  by  I.  Hence,  replacing  m  by 
I9+K-1,  we  arrive  at.  the  subscript 

KT  =  T  +  (I9+K-2) *MT 

for  T.  Here,  KT  is  the  field  subscript  which  refers  to  the  testing 
function.  By  analog> ,  the  source  subscript  LT  which  refers  tc  the  expan¬ 


sion  function  is  given  by 


LT  *  J  +  (.19 FI, -2)  *MT 


Retaining  as  the  testing  function,  we  take  the  analogously 

subscripted  function  fj^+^_^(t)  to  be  the  expansion  function  and  arrive 
at  the  matrix  subscript 

J1  =  (J9-E1-2)  *N2  +  I9+K-1 

where,  as  in  the  program,  N2  =  2*N, 

DO  loop  17  accumulates  in  R1  the 
equation  (31)  of  [91  due  to  the  - 

This  contribution  is  given  by 

2*MT 

R1  =  ^  1I+2*MT*(J-1)  *  TI+.-.*MT*(J-D*  pd(i+mt*(j~1>) 
where,  as  in  the  program, 


contribution  to  (Yt£)  of 
n  J  j  J 

term  in  equation  (32)  of  [9]. 


™(i)  ■  - H  *  Ero-) 

xi 

Here,  the  factor  -M  not  included  in  equation  (32)  of  [9]  provides  for 
the  choice  of  field  evaluation  either  outside  or  inside  S. 


DO  loop  18  accumulates  in  ..1  t  \e  conttibution  to  (Ytt:) 


O* 


n  'j-l,j 

equation  (Ji)  of  [9]  due  to  the  ,  -  tern  in  equation  (32)  of  19], 

■  ■  u .  p  . 

1  1 

This  contribution  is  given  by 


MT  ... 

R1  -  l  1  I  +2 *MT* (. f-  i, )  -MT  I  +2*M  1  I 

I  -1 


*  Pi)  ( I +KT*  (J-l ) ) . 


LISTING  OE  THE  SUBROUTINE  YZ 

SUBROUTINE  YZ  INN,  NIP,  N-*HI  »M,MTtRHtZH»X,a»Y,Z) 

COMPLEX  U»Y  <784) »Z  <  784  ) , Gi,G2tG3 ,G4,G5,G6, Y1 ,Y2, Y3,Y4,Zl ,Z2 , Z3 ,Z4 
PI  ME  NS  ION  KH 1161 ) , ZH ( 1 61 ) , X( 48 ) »  At  48 ) , D( 160)  , PD<  160) ,  TP (  32  0) 

PI  MENS  ION  CR  (4  8)  ,Cl(48  )  ,C2(4S)  ,C3(  48)  ,C4<48) 

COMMON  RSI  160)  ,ZS(  16  0)  ,  S  Y(  I  6  0)  » C  VI 1 6  0 )  ,  T  (  3  20  ) 

PI =3 .141593 

PI  m~~M*PI 

N=  (NP-1  )  /MT-1 

N2«2*N 

NG= I N+ I ) *MT 

NGM=NG~MT 

MT  2- MT  *2 

DC  5  7  1=1  r  NG 


12=1+1 

DR=RH(  12  )-*RH ( 1 ) 

DZ=ZH< 12 )— ZH( I ) 

D(  I » =SQRT(DR*OR+DZ*DZ> 

RSU  )=*5*(RH(  I2)+RH(I)  ) 
ZS(I  )  =  ,  5*(  ZH(  I  2)+ZH(  I)  ) 
SVU  )  =  DR/D (  I  ! 

CVU  )  =DZ/D ( 1 ) 

PPl I  )=  P IM/ ( D I  I  ) *RS (I ) ) 

5?  CONTINUE 
Jl  =  0 
J5=0 

00  6  8  J  =  1  » N 
0EL=0. 

DO  12  1=1, MT 
J1=J  1+1 
DEL=  DEL  +  DI J1  ) 

12  CONTINUE 
J 1  = J 1-MT 
SN=0  . 

00  19  1=1 ,MT 
J5  =  J5+  l 
J1  =  J  1+  1 
SN=SN+D( Jl) 

TPC J5)=D( Jl)  /DEL 
T(  J55=< SN-. 5*D< Jl)  )  *TP  ( J  5) 
19  CON V INU E 
DEL= 0. 

DO  15  I =1 ,MT 

J  1=  11*1 

DEL  =  DFL  +  D(  Jl  i 

15  CONTINUE 
J1  =  J  1-MT 
SN=DEL 

DO  16  I =  1 »  M  T 
J5-J6+  1 
Jl-J  1+  1 
$N= SN-Dt Jl  ) 

TP ( J6) =—  D(  Jl  1/OEL 
1  (  J  5  )  -  -  (  SN+  .  * D  (  J 1  )  )  *TP<  J5) 

16  CONTINUE 
Jl =J1-MT 

6  8  CC  NT  IN U I 
P 1  2  =  . ,  *  P 1 
EN-N  N 

DO  2  5  K "  1 , N P H I 
PH-P  i  2M  XI  K  >  »  1  .  ) 

PHN= PP# ( N 
SN--5  1  i.T'f'hl 
(  R  (  K  )  -  4  .  <  6  N  *  SN 
r  }  ,■  i  :•  *  a  i 1  ) 

C.S  --  E  1*C0S(  PHN) 

(  1  (  K  .  5*<  P  (K  )  *C  5 

<  .>  { K  }  <  ( ;  PP  \  *  (  .i 

('..?(  K  >  =  kl*SIN'<  I  N  (  P  HN) 


i iaaa&L 


C4<K  J  =  CS 
2  5  CONTINUE 
N2N=N2*N 
N4N=  N2N*2 
DO  62  J  -l  ,N4N 
Y(  J)=0. 

Z  i  J )  =0  . 

62  CONTINUE 
U=(0.,1  .  ) 

DO  59  J  =  1 »  NG 
FJ  =  FN/RS(J) 

L 1  =  1 
L2=2 

IFtJ.LE.MT}  L  1=2 
IF ( J  .GT . NGM }  L2=l 
J9=( J-l) /MT 
JT=J+MT*( J9-2) 

J5=( J9-2)*N2-1 
S 1  =  1 . 

DO  60  I=1,J 
I9=(  I-1I/MT 
I T  =  I  +M  T  *  (  I  9-2  i 
J6= I  9+ J  5 
FUFN/RSCI ) 

RP=RS( J)-RS(I ) 

ZP=ZS( J) -Z  S  C  I) 

R2=RP*RP+ZP*ZP 
IF (I .NE.J)  GO  TO  41 
SI  =.  5 

R  2= .  0625*0 ( J ) *D( J) 

41  R3=RS( I i*RS (J) 

GUO. 

G2=0. 

G3=0  . 

G4*0. 

G  5  -  0  . 

G  6  =  0  . 

00  6  1  K=  1  *  NPHI 
P4=R2+R3*CR  (K) 

R  5=  S  QR  T  (  R4  ) 

Z  1=  S )  /R  5*  (  C OS  (  R5 )-'J*  Si  N(R5)  ) 

Y 1=Z  1*( 1,  +  U*R5)/R4 
C  1  =  C  1  (  K  J  *Y  1  +  G1 
G2=C2( K 1 *V1 +G2 
G3 -C 3{ K ) *  V 1 +  G  3 
G4  =  C  4 1 K ) *Z  I+G4 
G5=C?(K  )  *  Z  1  +G5 
G6=C3(  K)  *Z  HG6 
61  CONTINUE 
GT=U*G3 

Yl  =  (  RP*C  V(  J  >-ZP*SV(  J  i )  *G2-RS  ( I  )  *CV  ».  J  )  *G1 

Y2=(  RS  <  J  )*S  V(  I  )*CV(  J1--RS  <  I  )*SV(Ji*CV(I  >-ZP*SV(  I5*SV(J))*G3 
Y  3=  Z P* G3 

Y4=  <  RP*CV(  I  )-2P*SV  (I) )  *G2+RS<  J  )*CVm*GI 
Z1  =  U*<  SV( I  U  SV< J)*G5+C VI 1 ) *CV ' J ) *G4 ) 

G1=-U*G4 
Z?  =  - SV( J ) *  G  6 
G2=~ FI*G4 
Z3=  S V ( I  } *G6 
G3=F  J*G4 
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Z4=U*<  G5--F  I*G3> 

Kl  =  1 
K2  =  2 

IFU.LE.MT)  Kl“2 
I F ( I »GT  .NGH }  K2  =  l 
00  3 1  L=L1 ,L2 
LT-JT+MT*1 
J7=J6+  L*N2 
00  32  K=K1,K2 
KT=I  T+MT*K 
TT=  T  ( L  T  ‘  ‘T(KT) 

J 1 =J  7+K 
J2=J 1+N 
J3=J  1+N2N 
J4=  J  3+N 

Y( Jl )=TT*Y1+Y< Jl) 

Y(  J2  f-^TT*Y2+Y<  J2> 

Y ( J3 )  =  TT :  Y3+Y( J3  » 

Yi J4)=TT*Y4+Y< J4) 

Z(  Jl  )=TT*Zl  +  TPaT)  *TP(  KT}*G1  +  Z(  Jl) 
Z< J2 )=TT*Z2+TP( LT)*T (KT) *G2+ Z ( J2 ) 
2<  J3)=TT*Z3+TP(KT)  *T  (L  T )  *C,3  +Z  i  J 3  S 
Z(  J4  )  =  TT*Z4*  Z( J4) 

32  CONTINUE 
31  CONTINUE 
60  CONTINUE 
59  CONTINUE 
N2P-N2+1 
KOI  =  1 
J1=0 
J5=0 

00  11  J=1,N 

R1  =  0. 

DO  17  I  =  1  -  MT2 

Jl=Ji+l 

J5=  J 5+1 

R1=R1+T(J1)  *T(  Jl  )*  PD  ( J  51 
17  CONTINUE 
J1=J 1-MT2 
J 5= J  5-MT2 
KD2  =  KO 1 +N 
K  03=  KOI  +N2N 
KD4= K03+N 
G1=Y (KOI )-Y (KD4) 

Y(KDl) -R l+Gl 
Y( KD2 )  =  0. 

Y  (  KD3 )  =  0  . 

Y<  KD4)  =R1— G 1 
Z<KD1)=Z<KD1)+Z<KD1) 

ZIKD2)  =Z(K02)-Z(  Kt)3J 
Z  <  KD3) =  -Z.  I KD2 ) 

Z(K04»=Z(K04»+Z(KD4) 

1 F ( J - I i  26,27,26 
27  J 1  =  J  1  +  M T 2 
J  5  =  J  5+  M  T 
GO  TO  22 
26  KU1=K01— 1 
KU2  =  K02— 1 
KU3=  KD3- 1 


,  'Mr*' 
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KU4*KD4-1 
KLl  =  KDi-N2 
KL2=  KD2--N2 
KL3  =  KD3-"  N2 
KL4=  KD  4— M2 
R1  =0  . 

00  18  1*1. MT 
J1=J1+1 
J5=J  5+1 

R1=R1+7(J1)*TI  Jl-MT)*FDt  J5) 
18  CONTINUE 
J1=J  1+MT 
G1=Y (KU1 S-Y( KL4) 

G2=Y (KU4  )-Y  ( KL  1 J 
Y( KU1) =R1+G1 

Y  (  K  U  2 ) =  Y ( KU2  5 - Y I KL2  5 
Y(  KU3  5  =  Y  ( K  U3  5  -*Y(  KL  35 
Y ( KU  4  5  =  R 1  +  G2 

Y  i  KL  1 5  =  R  1-G2 
Y( KL  2 ) =- Y( K  L2 5 
Y( KL  3) =-Y( KU3 ) 

Y  (  KL  4  5  =E  1-  G 1 
Z(KU1) =Z(KU1 )+Z(KLi) 
Z(KU2)=Z(KU2)-Z(KL35 
Z( KU3) =  Z ( KU3  5  ~Z ( KL2 ) 

Z(KU4) =  2 ( KU  4  5  +Z ( KL4 ) 

Z ( KL 1 )=Z (KU1 5 

Z( KL  2) =— Z( KU3 } 

Z ( KL  3  5  =-Z( KU2  5 
Z(KL4)  =Z ( KU4 ) 

22  KD 1  =  KD  1  +  N2P 
11  CONTINUE 

IF(N.LT<3)  RETURN 
J2  =  N2 

DO  13  1=3,  N 
J2  =  J  2+  N 2 

Jl=I-2 
KL  1=  I 

DO  14  J=1,J1 
KU1= J2+J 
KU2-KU1+N 
KU3=KU  1  +  N2N 
KU4=  KU3  +  N 
KL2- KL  1+N 
KL  3=  KL  1  +  N2N 
KL4--  KL.3  +  N 

Y  ( KL  ] )=-Y(KU4) 

Y  (  KL  2 ) =— Y { K U2 i 
Y( KL 3 5  =— Y (  KL3  5 
Y( KL  4 ) =-YC KU1) 

Z(KL 1I-ZIKU1! 

Z  ( KL  2)  =—Z  (  KU  3  ) 
Z<KL3)=-ZtKU2» 

Z  (  KL  4 )  =  Z  ( K  U4  J 
KL 1= KL 1+N2 
14  CONTINUE 
13  CONTINUE 
RETURN 

end 
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III.  THE  SUBROUTINE  PLANE 


Description: 

The  subroutine  PLANE (NN,  N,  MT,  NT,  THR,  R)  puts  of  (49)  In 
R(J  +  (w-l)*N+(L-.l)*4*N)  where  j*=l,2,,..N,  m=l  denotes  pq  =  t9,  m**.: 
denotes  pq  »  <|>e*  m=3  denotes  pq  =  t<J>,  m=4  denotes  pq  -  <jx(>,  and  L«*l,2, .  .  .NT. 
Here,  L  denotes  the  Lth  value  of  the  receiver  angle  9  .  The  first  5 
arguments  of  PLANE  are  Input  variables.  Except  for  the  new  input  v*  ‘table 
MT,  the  subroutine  PLANE  is  the  same  as  the  old  subroutine  PLANE  on  pages 
22-26  of  [10].  If  MT  =  2,  these  subroutines  are  exactly  the  same  as  Car 
as  the  calculation  of  R  in  terms  of  the  rest  of  the  input  variables  is 
concerned. 


The  integration  of  Ct)  over  t  inherent  in  (49)  is  approximated 

by  satnp’ing  f_.  (t)  2*MT  times  instead  of  4  times.  The  representation  of 

pfc^Ct)  given  by  (66)  of  [9]  is  replaced  by  (77).  NN  is  the  value  of  n 

appearing  in  R^.  It  is  required  that  NN  >_  0  but  this  requirement  causes 

no  real  loss  of  generality  because  rP^  is  either  even  or  odd  in  n,  N  is 

n3 

the  number  of  expansion  functions  lying  cn  the  generating  curve.  Specifi¬ 
cally,  N  is  the  maximum  value  of  i.  in  (77).  THR(L)  is  the  Lth  value  of 
the  receiver  angle  where  L  =  1,2,... NT.  The  variables  RS,  ZS,  SV,  CV, 
and  T  appearing  in  the  common  statement  e^rly  in  the  subroutine  PLANE  are 
input  variables  calculated  by  calling  the  subroutine  YZ  beforehand.  The 
calculated  values  of  these  variables  depend  only  on  the  second,  fifth, 
sixth,  and  seventh  arguments  (NP,  MT,  RH,  ZH)  of  YZ. 

Minimum  allocations  are  given  by 

COMPLEX  R(4*NT*N) 

DIMENSION  THR(NT) ,  BJ(M) 

COMMON  RS(NG) ,  ZS (NO) ,  SV(NG),  CV(NG) ,  T(2*MT*N) 

where 

NG  =  (N+1.)  *M  i 


and  M  is  the  largest  of  the  values  of  M  calculated  by  PLANE.  The  sug¬ 
gested  allocation  Bd(fO)  will  work  if  the  maximum  circumference  of  the 
body  of  revolution  is  less  than  26  wavelengths. 
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Most  of  the  statements  in  the  subroutine  PLANE  are  the  same  as 
or  very  similar  to  statements  in  the  old  subroutine  FI.ANE  listed  on 
pages  25-26  of  [10].  The  major  difference  is  in  the  calculation  of 
the  subscript  (TT+MT*K)  for  T  and  the  subscript  J1  for  R,  Using 
reasoning  similar  to  that  used  to  obtain  the  subscripts  KT  and  J1  in  the 
subroutine  YZ,  we  arrive  at 


(IT  +  MT*K)  =  I  +  (19  +  K-?)*MT 
J1  =  19  +  K-l 

where 

19  =  [(I-D/MT] 

The  above  Jl  is  valid  only  for  L=l.  If  L  ^  1,  then  (L-1)*4*N  must 
be  added  to  this  Jl. 


LISTING  OF  THE  SUBROUTINE  PLANE 

S’JBROU'  INF  <>LANE(NN,  N,  MT,NTt  THR  »  R ) 

COMPLEX  R1  1064)  »  U,  U1  .0  2*  Ri  *R2  ,R3  *R4 
DIMENSION  THRI37) ,  BJ  (50) 

COMMON  RS(  160)  »ZS(  160)  ,SV(  160)  ,C  V(  160)  ,T  (320) 
NG=( N+l )*MT 
U=  (  0  .  t  1  .  ) 

Ul=3 ,141E93*U**NN 

N4=4*N 

Jf’~N4*NT 

D(  2?  J=1,JF. 

R  (  J  ')  -0. 

22  continue 

J  5  =  -  1 

DO  12  =  1  *  N  i 

cs-ccs  (Thp  im 

SN-2 •“SI N ( THM(L) ) 

DO  13  1=1 .MG 

X  -  .  ?  5  *-  k  S  (  I  )  *  S  N 

IHX.t.  F..5E  7)  Gil  TO  1C 

M-  2 .  8*X  +  13.-2-  /x 

IF(  X.LT..5)  m=10.6  :  AF  "G10(X) 

I  F  (  M  .(» F  .  (  NN  +  2  )  )  GO  10  19 

18  B  J I  -  0 . 

BJ2  =  0. 

B  J 3 -  0. 

if(Nn.;-:c.i  )  b  j i  =  i . 

IF(NN-EO.O)  B.i/'-l. 

go  rn  24 

19  B  J  t  M  )  -  0  <• 

J  M  =  M  - 1 

Bj  < J  M } -  i  • 
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00  14  J - 3  t  M 
JM=  J  M— 1 

BJ(  JM)  - JM/X*B  J (  JM+ 1 J  ~S  J (  JM+2  ) 

14  CONTINUE 

s^o« 

DO  15  J  - •  3  ,  M » 2 
S=S+BJ(  J) 

15  CONTINUE 

S~ BJ  ( 1 )+2.*S 
BJ2=BJ  (NN+1) /S 
8  J3~  BJ ( NN  +  2 I / S 
B J1 “-B  J3 

IF(NN.GT.O)  BJ1=BJ(NN)/S 
24  ARG=ZSUJ*CS 

U2=U1*(CCS< ARG)+U*SIN< ARG) i 
R4=<  BJ3-BJ1  »*U*U2 
R2  =  (  BJ 3+B J 1 ) *U2 

R 1=~  BJ2*CV ( I )*SN*U2  +  CS#SV( I ) *R4 

R3=S  V ( I  )  *R2 

R2=~CS*P2 

1 9=  (  1-1)  /MT 

I T= I +MT  1 9-2  I 

J7= I  9+ J  5 

K1  =  1 

K2=2 

If(I9»EQ«0)  KI-? 

I F ( J9.EQ.N)  K  2 - 1 
DC  20  K=K1,K2 
TT-’T  (IT+MT*K) 

Jl=  J7+K 
J2  =  Jl+N 
J3  =  J  2+N 
J4= J  3+N 

R ( J 1 ) =TT*Ri+R{ Jl) 

R(  J2  )■-■  T T*R  2*R  I  J2 ) 

R  ( J3 )»TT*R3+R( J3) 

R ( J4 ) =  T  T*R4  +  R ( J4 ) 

20  CONTINUE 
13  CONTINUE 
J 5= J  5+  N4 
L  2  CONTINUE 
RETURN 
END 
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IV.  THE  SUBROUTINES  OECOMP  AND  SOLVE 


Description : 

The  subroutines  DECOMP (N, IPS, UL)  and  SOLVE (N, IPS, UL,B,X)  solve 
a  system  of  N  linear  equations  in  N  unknowns.  These  subroutines  will 
be  used  in  Section  V  to  solve  the  matrix  equation  (39).  The  input  to 
DECOMP  consists  of  N  and  the  N  by  N  matrix  of  coefficients  on  the  left- 
hand  side  of  the.  matrix  equation  stored  by  columns  in  UL.  The  output 
from  DECOMP  is  IPS  and  UL.  This  output  is  fed  into  SOLVE.  The  rest  of 
the  input  to  SOLVE  consists  of  N  and  the  column  of  coefficients  on  the 
right-hand  side  of  the  matrix  equation  stored  in  B.  SOLVE  puts  the 
solution  to  the  matrix  equation  in  X. 

Minimum  allocations  are  given  by 

COMPLEX  UL(N*N) 

DIMENSION  SCL(N) ,  IPS(N) 

in  DECOMP  and  by 

COMPLEX  UL(N*N) ,  B(N) ,  X(N) 

DIMENSION  IPS (N) 

In  SOLVE. 


More  detail  concerning  DECOMP  and  SOLVE  is  on  pages  46-49  of  [13] 

LISTING  OF  THE  SUBROUTINES  DECOMP  AND  SOLVE 

SUBROUTINE  DECPMPC N* IPStUL) 

COMPLEX  UL (3136) »P1 VOT.EM 
DIMENSION  SCL(  56) ,  I  PS(  56) 

DO  5  I  -  I » N 
I  PS ( I )=  I 
RN  =  0 . 

JI=I 

DO  2  J  =  1 1 N 

ULM= ABSjRE AL ( UL(JI ) )  )  +  AB $C  A I  MAG < UL  t J 1 ) l ) 

J1=J 1+N 

IF(RN-ULM)  1  »2  »2 

1  RN=ULM 

2  CONTINUE 
SCL ( 1)=1 ./RN 

5  CONTINUE 
NM 1  =  N-  1 
K2--0 

DO  IT  K= 1»  NM1 
B I  G=  0. 

DO  11  I  *  K  ,  N 
I  P=  I  PS  U  i 
1PK= IP+K2 
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SIZE=C  ABS<  REAL  (ULl  EPK1  ) )  ♦ABS  t  AIMAG  I  UL  (  IPK )  ) )  )*SCU1P) 
IF(SIZE-BIG)  ll»U*iO 

10  B I  G=  S I Z  E 
IPV=  2 

11  continue 

IFI  I  PV--K3  1 4 » 1 5 *  J. 4 
14  J=  I P  SC  K ) 

IPS! K)  =  IPS ( IPV ) 

I  PS  C  1 PV ) =J 
1  5  KPP=IPS<  K)  +  K2 
P I  VO  T~Ul.  (  KPP  ) 

KPl-K+l 

DO  16  I =KP1 f N 

KP-KPP 

I P  =  I  P S  <  I  )  +  K2 
EM~-UL<IP)/PIVOT 
18  UL  ( I  PI  =— EM 

DO  16  J~  KP 1 »  N 
IP- I P+N 
KP-K  P+N 

ULU  P)-UL<  IP)  +  EM*UL<  KP) 

16  CONTINUE 
K2*K2+N 

17  CONT  INUE 
RETURN 
END 

SU6R0UT INE  SOLVE ( N, I  PS ,Ul ♦ 8 » X ) 

COMPLEX  UL  131365  ,B(  56)  ,X( 56  )  ,  SUM 
DIMENSION  IPS!  56) 

NP1  =  M+  1 
IP  =  IPS(  1  i 
X( 1 ) =B( IP) 

DO  2  1=2, N 
IP=I PS( I) 

IPB  =  IP 
1  M 1  =  1  - 1 
SUM  =  0. 

DC  1  J  =  1  , IMi 

SUM  =  SUM-s  UL  (  IP)  ♦XIJ) 

1  IP=IP*N 

2  XI  I  )  =B  I IPB  »~SUM 
K2=N*(  N~  1) 

IP- I  PS  tM)  +  K2 
XI  N)  =X(  N.i/ULJ  I  P) 

DO  &  I  BACK*?,  N 
I=NP1-I BACK 
K2=K2-N 
I P  I  =  I P  S I  I )  +  K  2 
IP  1=  1*1 
SUM*  0. 

I  P=I  PI 

DO  3  J  =  I P 1  iM 
IP= I P+N 

3  SUM=  SUM  .  UL.  ||P|*XIJ) 

4  XU  i  -  (  X(  D*  SUM  )/UL  I  I  PI) 

RETURN 

END 
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V.  THE  MAIN  PROGRAM 


Description: 

The  main  program  calculates  the  electric  and  magnetic  currents 

^00  ^i0 

(58)-(59)  and  the  normalized  scattering  cross  sections  jy-  and  -*7  (70) 

for  the  0  polarized  axially  incident  (8  =  180°)  plane  wave  (50) — (5 1) . 

The  main  program  calls  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE. 

The  main  program  is  short  and  simple.  It  is  a  representative  applica¬ 
tion  of  the  theory  in  Part  One  of  this  report. 

Input  data  is  read  early  in  the  main  program  according  to 

READ(1,10)  NP,  NPHI ,  MT,  NT 

10  FORMAT (413) 

READ (1,11)  BK,  UR,  ER,  ALP,  BET 

11  FORMAT (5E14. 7) 

READ(1,12) (RH(I) ,  1=1,  NP) 

READ(1, 12) (ZH(I) ,  1=1,  NP) 

12  FORMAT ( 10F  8 . 4) 

READ (1, 11) (X(K) ,  K=l,  NPHI) 

READ(1, 11) (A(K) ,  K=l,  NPHI) 

The  input  variables  NP,  NPHI,  MT,  RH,  ZH,  X,  and  A  are  very 
similar  to  variables  of  the  same  names  in  the  argument  list  of  the  sub¬ 
routine  YZ.  In  summary,  (RH(i) ,  ZH(i)),  1=1, 2,... NP,  are  the  cylindrical 

coordinates  (p^,  z^)  on  the  generating  curve,  2*MT  is  the  number  of  values 
of  t  at  which  f .  (t)  is  sampled  for  the  purpose  of  numerical  integration, 
and  X  and  A  are  respectively  the  abscissas  and  weights  for  the  NPHI  point 
Gaussian  quadrature  integration  in  <J>. 

The  scattering  cross  sections  are  evaluated  at  receiver  angles 

0r  =  (J-lW(NT-l)  radians  for  J=1,2,...NT.  BK  is  the  propagation 

constant  k  in  the  external  medium.  This  k  appears  in  (50)-(51).  UR 

pd 

and  ER  are  respectively  the  relative  permeability  —  and  the  relative 

j 

permittivity  —  of  the  body  of  revolution.  Here,  and  are  the 
£e 

permeability  and  permittivity  of  the  (diffracting)  body  of  revolution 
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and  u  and  e  are  those  of  the  external  medium.  ALP  and  BET  are  respec- 
e  e 

tively  the  constants  a  and  6  appearing  in  (39).  The  PMCHW  solution  is 
obtained  by  setting  a  =  8  =  1.  The  Muller  solution  is  obtained  if  a  and 
$  are  given  by  (22)  and  (23) . 

Minimum  allocations  are  given  by 

COMPLEX  YE (4*N*N) ,  ZE(4*N*N),  R(4*NT*N) , 

B(4*N),  YD(4*N*N) ,  ZD(4*N*N) „  Y(16*N*N), 

C(4*N) 

DIMENSION  RH(NP),  ZH(NP),  X(NPHI),  A(NPHI), 

THR(NT),  RC(N) ,  IPS(4*N) 

where 

N  =  [(NP-1)/MT]  -  1 


Statement  38  puts  Y^e  defined  in  accordance  with  (75)  by 


,tt 

Yt+ 

le 

le 

•<frt 

le 

le 

in  YE  and  Z defined  in  accordance  with  (76)  by 


(84) 


Z 


le 


(85) 


in  ZE.  Storage  of  Y^e  and  Z^g  is  by  columns. 


Statement  39  puts  the  matrix  R^  defined  by 


Rl  = 


*te 


*40 

R1 


*1 


(86) 


in  R.  The  column  vectors  on  the  right-hand  side  of  (86)  appear  in  (56) . 
For  receiver  angle  6^  =  THR(J),  R^  of  (86)  is  stored  by  columns  in 
R((J-1)*4*N+1)  through  R((J*4*N). 
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DO  loop  22  uses  (56)  to  store  the  right-hand  side  of  (39)  in  B. 

DO  loop  22  also  puts  2/(kp_ . )  in  RC  where  pTXl_. ,  is  the  cylindrical 
coordinate  radius  evaluated  at  the  peak  of  the  triangle  function  inherent 
in  fI(t) . 

Statement  40  puts  defined  by 


(87) 


in  YD  and  defined  by 


(88) 


in  ZD.  Storage  of  Y^  and  Z^  is  by  columns. 

Nested  DO  loops  26  and  27  put  the  first  two  columns  of  submatrices 
on  the  left-hand  side  of  (39)  in  Y.  The  index  J  of  DO  loop  26  denotes 
the  Jth  column  of  the  composite  square  matrix  on  the  left-hand  side  of 
(39). 


Nested  DO  loops  28  and  29  put  the  third  and  fourth  columns  of  sub¬ 
matrices  on  the  left-hand  side  of  (39)  in  Y.  The  index  J  of  DO  loop  28 
denotes  the  (2*N+J)th  column  of  the  composite  square  matrix  on  the  left- 
hand  side  of  (39). 


Statements  41  and  42  solve  the  matrix  equation  (39)  for  the  com¬ 
posite  column  vector  consisting  of  V^,  V^,  1^,  and  1^.  This  composite 
column  vector  is  stored  in  C. 

At  the  peak  of  the  Jth  triangle  function,  the  n=l  term  of  the 
equivalent  electric  current  (58)  reduces  to 
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i  “  ^^pmt*j+i^  xij  Et008*'1'  ^2^pmt*j+i^  zij  ^>sin^  <'89^ 

DO  loop  31  prints  the  real  and  imaginary  parts  of 

U1  =  (2/(kPMT*j+1))I1j  » 
the  real  and  imaginary  parts  of 


U2  =  (2j/(kp‘T*J+1»I  l 

and  the  magnitudes  of  U1  and  U2.  Here,  U1  is  the  t  component  of  the 
equivalent  electric  current  in  the  E  plane  and  U2  is  the  4>  component 
of  the  equivalent  electric  current  in  the  H  plane  when  the  y  component 
of  the  incident  magnetic  field  is  minus  one  at  the  origin. 

At  the  peak  of  the  Jth  triangle  function,  the  n=l  term  of  the 
equivalent  magnetic  current  (59)  reduces  to 

£  /  Wi> vu  V1"  * +  <2/pmi*j+i'  vu  V  °s  *  < 

DO  loop  34  prints  the  real  and  imaginary  parts  of 


BCJ)  -  «J/0»W.m»V1.7 
and  the  real  and  imaginary  parts  of 


“  '  <2/<»Wl))VU 

and  the  magnitudes  of  B(J)  and  U2.  Here,  B(J)  is  the  t  component  of 
the  equivalent  magnetic  current  in  the  H  plane  and  U2  is  the  <f>  component 
of  the  equivalent  magnetic  current  in  the  E  plane  when  the  x  component 
of  the  incident  electric  field  is  minus  one  at  the  origin. 

Nested  DO  loops  35  and  36  calculate  and  print  -p-  and  -p~  °f  (70). 
Inner  DO  loop  36  accumulates  the  portion 
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+  R**S*e  +  rJ8^8 


+ 


of  (62)  in  ET  and  the  portion 


-i^;8  -r*8?;8  +  r[*?;9  +  sj*^8 


of  (63)  in  EP.  The  W1  and  W2  printed  in  DO  loop  35  are  respectively 


aee 


^  and  -^r-  for  receiver  angle  0  =  (J-l)*ir/(NT-l)  radians. 


Suggested  modifications  of  the  main  program  are: 


1)  Changing  the  normalization  of  the  scattering  patterns. 

For  example,  one  could  replace  (70)  by  (69).  All  scat¬ 
tering  patterns  in  Part  One,  Section  V  are  plots  of  (69). 

2)  Removing  the  restriction  that  the  values  of  the  input 
arguments  NP,  NPHI,  MT,  RH,  ZH,  X,  and  A  of  the  subroutine 
YZ  be  the  same  in  call  statements  38  and  40.  This  modifi¬ 
cation  is  indicated  by  (74). 

3)  Generalizing  from  axial  plane  wave  incidence  to  oblique 
plane  wave  incidence. 


The  above  three  modifications  can  be  realized  without  tampering  with 
any  of  the  subroutines  YZ,  PLANE,  DECOMP,  and  SOLVE. 

The  sample  input  and  output  accompanying  the  listing  of  the 
main  program  is  for  the  dielectric  sphere  with  ka  =  1  and  =  4. 
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LISTING  OF  THE  MAIN  PROGRAM 

//PGM  JOB  <XXXX,XXXX,1,2) ,*MAUTZ,JOE* ,PEGI0N=200K 
//  EXEC  WATFIV 
//GO  .SYS  IN  OD  * 

SJOB  M AUTZ , TIME= It  PAGES=40 

C  SUBROUTINES  YZ,  PLANEt  OECOMP,  AND  SOLVE  ARE  CALLED. 

COMPLEX  YE  (784)  ,ZE  (784  ),  R  (1064)  ,0(56)  ,  YD(784  ) ,  ZD (  784  ) ,  Y ( 3 136  ) 
COMPLEX  C(56),U,UltU2,ET,EP, CONJG 

DIME  NS  ION  RH(161) , ZH<161 ), X(48) , A(  48) ,THR(37 ) ,RC ( 19 ) . I  PS ( 56) 
READ(ltlO)  NP , NPHI , MT, NT 

10  FORM  AT  (413) 

READ (1  ,11)  BK,UR,ER,ALP,  BET 

11  FORMAT ( 5E14.7) 

READ(1,12) (RH( I) ,1=1, NP) 

READ ( 1 , 12 ) ( ZH( I ) , I  =  1 ,N  P) 

12  FORMAT ( 1  OF 8.4 ) 

READ(1,11)(X(K),K=1,  NPHI  ) 

READ(1,  11)(  A(K),K=1,NPHI  ) 

WRIT  E( 3 , 13 )  NP,NPHI,MTtNT 

13  FORMATC  NP  NPHI  MT  NT*  /( 1 X  ,  13  ,  15 ,2  13  )  ) 

WRIT  E( 3 , 14 )  BK,UR,ER,ALP,BET 

14  FORMAT  (7X,'BK*  ,  12X  ,  *  UR  '  ,  12X  , •  ER  •  ,  1 IX,  •  AL  P*  ,1 1 X,  *  BET*  /  ( 1  X,5  E  14.7 )  ) 
WP I TE( 3 ,  15) (RH(I ) ,1  =  1, NP) 

15  FORM  AT  (  *  RH'/(  1X.10F8.4)  ) 

WRI TE( 3 , 16) ( ZH( I ) , 1  =  1, NP) 

16  FORMATC  ZH ' /(  IX ,  1  OF  8.  4)  ) 

WRITE(3,17)(X(K) ,K=1,NPHI) 

17  FORMATC  X* /(  1 X , 5E14 .7  ) ) 

WRIT E(  3,18) (A(K) ,K=1,NPHI) 

18  FORMATC  A* /(  1X,5E14.7  )  ) 

DO  19  J=  1, NP 
RH(J)=BK*RH(J) 

ZH( J  )  =  BK*ZH ( J) 

19  CONTINUE 

38  CALL  YZ(1,NP,NPHI,-1,MT,RH,ZH,X,A,YE,ZE) 

WRITE (3,20) (YE( I), 1=1,4) ,(ZE(I) , 1=1,4) 

20  FORMATC  SOME  ELEMENTS  OF  YE  AND  ZE  •  /  ( IX,  4E14. 7)  ) 

PI=3. 141593 

OT=P  1/ (NT— 1 ) 

00  21  J  =  1 ,  NT 
THM  J)«(  J-1)*0T 

21  CONTINUE 

N=( NP- 1 ) /MT— 1 

39  CALL  PLANE( 1 , N ,MT, NT ,THR ,R ) 

N2=2*N 

N3=3*N 

N4=4*N 

NTN= ( NT- 1)  *N4 
00  22  1=1, N 
J3=I+NTN 
B  ( I )  =R  (  J  3 ) 

B( I+N) =— R ( J  3+N ) 

B ( I+N2 )=-R< J3+N3) 

B( I+N3)=-R( J3+N2) 

RC ( I )=2./RH(MT*I+l) 

2  2  CONTINUE 

WRITE! 3,23) (B( I) ,1=1 ,N2) 

23  FORMATC  HALF  GF  THE  ELEMENTS  OF  B'  /  ( 1 X  ,4E  14 .7 )  ) 

EM=SQRT ( UR*ER) 
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DO  24  1=1, NP 
KH(I  )=EM*RH(I  ) 

ZH(I )=EM*2H( I) 

24  CONTINUE 

40  CALL  YZ(1,NP,NFHI,1,MT,RH*ZH,X,A,YD,ZD) 

HRlTE(  3,25)  (YD!I)*I=1,4)  ,  ( ZD ( I )  ,1=1,4) 

25  FORMAT ( »  SOME  ELEMENTS  OF  YD  AND  ZD'/! 1X,4E14.7)) 

D=  SQRT ( UR/ ER ) 

ALPD=ALP*D 

BETD=BET/D 

JY=0 

J1=0 

DO  26  J=1,N2 
DO  27  1=1, N 
JY=J Y+ 1 
J1=J1+1 
J2=J  1+N 

Y !  JY  )=YE(J2)+ALP*YD(J2J 
Y!  JY*N ) =— YE ( J1 }-ALP*YD (J 1) 

Y(  JY+N2 ) =ZE (J2)  +  BETD*Z0< J2) 

Y(JY+N3)=-ZE! J1)-BETD*  ZD(J1) 

27  CONTINUE 
JY=J Y+N3 
Jl=J  1+N 

26  CONTINUE 
J1=0 

DO  28  J=1 ,N2 
DO  20  1=1, N2 
JY= J Y+ 1 
Jl=Jl+l 

Y!JY)=ZEIJ11+ALP0*Z0(J1) 

Y(  J  Y+N2 ) =YE ( J1 )  +  BET*YD !  J1 ) 

29  CONTINUE 
JY=JY+N2 

28  CONTINUE 

41  CALL  DEC0MP(N4,IPS,Y) 

42  CALL  SOLVE! N4, IPS, Y,B,C) 

WRITE! 3,30) 

30  FORMAT! *0  REAL  JT  IMAG  JT  REAL  JP  IMAG  JP  MAG  JT 
1  MAG  JP') 

U=(0.,1.) 

DO  31  J  =  1 7  N 
RR=RC( J I 
ET=U*RR 
B( J) =  ET*C ( J ) 

J1=J+N 

B ( J 1  )  =  RR*C(J1) 

U1“RR*C ( J«-N2 l 
U2=E  T*C ! J+N3 ) 

W1=CABS(U1) 

W2=CABS<  U2  ) 

WP  IT  E(  3 ,32  )  U1  ,U2 , Wl  ,W2 

32  FORMAT! 1 X,6E11. 4) 

31  CONTINUE 
WRITE! 3, 33) 

33  FORMAT! *0  REAL  MT  IMAG  MT  PEAL  MP  IMAG  MP  MAG  MT 
IMAG  MP") 

DO  34  J  =1 , N 
U2  =  B  <  J  +N ) 

W1=CABS!B( J)  ) 
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W2=CABS(U2) 

WR  IT  E(  3  ,32 )  B(J)»U2,W1»W2 
34  CONTINUE 
WRIT  E(  3 , 37 ) 

3  7  FORM  AT  (  •  0  SIGTHETA  SIGPHIM 
CON= .25/PI**3 
JR  1  =  0 

DO  35  J  =  1 » NT 
ET=0. 

EP=0  . 

DO  36  1=1, N2 
JC2=  I+N2 
JR1= JR1+1 
JR2= JR 1+N2 

ET=ET+R(JR2 )*C(I)+R( JR1)*C( JC2> 

EP=EP-R(JR1  )*C(I)+R( JR2)*C(JC2) 

36  CONTINUE 
JR1= JR1 +N2 
W1=C  ON*ET*CONJG( ET ) 

W2=C  CN*E  P*CCNJG( EP ) 

WRIT  E( 3 , 32 )  W1,W2 
35  CONTINUE 
STOP 
END 

$DAT  A 

21  20  2  19 

0.1  000000 E+0 1  0.  ICOOOOQE+Ol  0.4000000E+01  0.  100Q000E+01  0.10000C0E+01 
0.  0000  0.1564  0.3090  0.4540  0.5878  0.  7071  0.8090  0.8910  0.9511  0.9877 

1.  0000  C.  9877  0.9511  0.  8910  0.8090  0.7071  0.5878  0.4540  0.3090  0.  1564 

0.  0000 

-1.0000  -0.9877  -0.9511  -0.  8910  -0.8090  -0.7071  -0.5878  -0.4540  -0.3090  -0.  1564 
0.0000  0.  1564  0.3090  0.4540  0.5878  0.7071  0.8090  0.8910  0.9511  0.9877 

1.  0000 

-0. 993 1286  E+ 00-0. 9635719E+ 00-0. 9122344F.+  00-0. 83911 70E+00-0. 74633 19  E+00 
-0.6360537E+00-0. 5 10 8670E+ 00-0. 3737061 E+00-0.2277859E* 00-0. 76526 52E-01 
0.7652652  E- 01  0 . 22778 59E+ 00  0.3737061E+00  0.  51 0 8670E+ CO  0.63605 37 E+00 
0.7463319E+00  0. 83911 70E+00  0. 912 2344E+00  0.9639719E+00  0.9931266E+00 
0 . 1 76 1401 E-0 1  0. 406 01 43 E- 01  0.6267205E-01  0. 8327674E- 01  0 . 1 0193  Cl E+00 
0.1 181945E+00  0. 1316886E+00  0. 1420961 E+00  0 . 149 17 30E+00  0 . 1527534E+00 
0. 1527534E+00  0.  1491730E+ 00  0.  142 096 1F.+  0 0  0. 1316886E+00  0 . 1 181945 E+00 
0. 1 019301 E+00  0. 8327674E-01  0 .626 720 5 E-01  0. 4060143E-01  0.1761401E-01 
$STOP 
/* 

// 

PRINTED  OUTPUT 
NP  NPHI  MT  NT 
21  20  2  19 

BK  UR  ER  ALP  BET 

0.1 090000E+  01  0. lOGOOOOE  +  Ol  0.4000000E  +  01  0 . 1000000E+01  0.1000000E+01 
R  H 

0.0000  0.1564  0.2090  0.4540  0.5878  0.  7071  0.809C  0.8910  0.9511  0.9877 

1.  0000  0.9877  0.9511  0.  8910  0.8090  0.  7071  0.5878  0.4540  0.3090  0.1564 

0.0000 
ZH 

-1.0000  -0.9877  -0.9511  -C.  8910  -0.8090  -0.  7071  -0.5878  -0.4540  -0.3090  -0.1564 
0.  0000  C.  1 564  0.3090  0.4540  0.5878  0.7071  0.8090  0.8910  0.9511  0.9877 

1.  0000 
X 

-  0. 993 1 286E+ 00- 0. 963 971 9E+ 00-0. 912 2344E+00-0. 8391 170E+00-0. 74633 19 E+00 
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-0.6 36 0537E+0 0-0.510 86 70E+  00-0.373  706 1 E+00-0 .  2277859E+00-0 .76526  52E-01 
0.7652652E-01  0. 2277859E+00  0.3737061E+00  0. 5108670E+ 00  0.6360537E+00 
0.7463319E+00  0. 83911 70E+ 00  0.9122344E+00  0.96397  19E+00  0  .993 1286E+00 
A 

0.1 76 1401 E-01  0.4060143E-01  0.6267208E-01  0. 8327675E- Cl  0. 1 0193C1 E+00 
0.1 181945F+00  0. 1316886E+C0  0. 1420961E+00  0. 1491730E+00  0. 1527534F+00 
0.1 52  7534E+00  0. 1491730E+00  0.1420961E+00  0. 1316886E+00  0. 1181945E+00 
0. 1019301 E+00  0.  8227675E-01  0.6267208 E-01  0.4060143E-01  0 . 176 14C1E-01 
SOME  ELEMENTS  OF  YE  AND  ZE 

0.2431891 E+01— 0. 5483869E—02  0 .5072693E+00-0. 1154096E-01 
0.  1 01 5252  E  + 00-0. 2 041  826 E-01  0.917  86  27E-01-0.30667  00E-01 
0.9008127 E- 01-0. 7025097 E+01  0.7607359E-01  0. 1940044E+01 
0.5597932E-01  0. 7892063E+00  0. 3345479 E-01  0 .2816698 E+00 
HALF  OF  THE  ELEMENTS  OF  B 

-0.5 464500 E+00-0. 7503574E+00-0.5449 371 E+00-0. 56906 88E+ 00 

-  0.4 7095 02 E+ 00-0. 3229641 E+ 00-0. 28 08 148 F+00-0. 1054077E+00 
0. 0 000000 E+00-0. 1778790E- 01  0.2808146E+00-0. 1054077E+00 
0.4 709504E+ 00-0. 3229641 E+00  0. 54493 70E+00-0. 5690686 E+00 
0.5 46 45 02 E+00— 0.7 50 35 76 E+00  0.7951685E+00-0. 5814634E+00 
0.7044444E+00— 0. 68  536  59E+00  0 .5376098E  +  00-0.  8194968E+00 
0. 2 93 905 8E+ 00-0. 93271 49E+ 00-0. 0000000 E+00-0. 976 96 01 E+00 

— 0.2 939057E+00-0. 93 271 47E+ 00- 0.53 76099 E+00-0. 8194970E+00 
-0. 7 04 4442 E+00-0. 68536 57E+00-0. 79 51688E+00-0.5814636E+00 
SOME  ELEMENTS  OF  YD  AND  ZO 

—  0.22  9096 9 E  +  01— 0 .37680 66E- 01 -0.3189305E+00—0.7057631E-01 
0.9 59 9513 E- 01-0. 1088477E+00  0.49421 2 5 E-0 1-0. 1374637E+00 
0.3 42 89 59 E+00— 0.29 021 17 E+01  0.2818084E+00  0. 1262741E+C1 
0.2  08  2631 E+00  0. 5027283E  +  00  0. 1465860E+00  0 . 19 1 1734E+00 


REAL  JT  IMAG  JT  REAL  JP 
-0.6398E+00-0.1395E+01  0.6170E  +  00 
-0.8349E+00-0.1227E+G1  0.8652E+00 
-0.1025E+01-0.9274E+00  0.728 5E +00 
-0.1127E+0 1-0. 4861 E+00  0.6632E+00 
—0.1 094E+0 1  0.5150E-01  0.4185E+00 
-0.8 731 E  +  00  0. 5874E+00  0.1349E  +  00 
—0.541 OE+OO  0.1026E+01  0.1116E+00 
-0.162 6E+00  0. 1320E+01-0.9319E-01 
0.1 722E  +  0 0  0. 1479E  +  01  0.2244E  +  00 


IMAG  JP 
0. 1415E+01 
0. 11 74E+01 
0.9655E+00 
0. 6917E+00 
0. 63 16E+00 
0. 71 85  E+00 
0. 1011E+01 
C. 1256E+01 
0. 1496E+01 


MAG  JT 
0.1 53  5E+01 
0. 1484E+01 
0. 1382E+0 1 
0. 1228E+0 1 
0.1 096E+01 
0.1052E+01 
0.  1 160E+01 
0.133  OE+0 1 
0. 1489E+01 


MAG  JP 
0.1543 E+Cl 
0. 1459  E+01 
0.1210E+Ci 
0.5582E+00 
0.7576E+00 
0.7311 E+00 
0.  1017E+01 
0.1260E+C1 
0. 1513E+01 


REAL  MT 
-0.1  064E  +  00-0 
•0.2323E+0  0-0 
-0.3712E+00-0 
-0.4731E  +  00-0 
-0.5  072E+00  0 
-0.4451E+00  0 
-0.3 171E+00  0 
-0. 1 556E+  00  0 
-0.1  304E-0  1  0 


IMAG  MT 

•  68C6E  +  00- 
.  5834E+C0- 

•  4190E+00- 
.  1809E  +  C0- 
.  1093E  +  00- 
.  399 2E  +  C0 
-6397E+00 
.  8002E+00 
.  8852E  +  C0- 


REAL  MP 
0.8434E-01- 
0.2176E  +  00- 
•0.163  9E+00- 
0.1442E  +  00- 
0.4617E-01- 
0.535  1E-01- 
0.701  8E-01- 
0 . 1 272E+00* 
0.1307E-01- 


I MAG  MP 
0.6963E+00 
0.5  837E+00 
0. 5138  E+00 
0. 42  09E  +  00 
0.  4199E  +  00 
0. 4940E+00 
0.6418E+00 
0.  7901  E+00 
0.  89  71E+00 


MAG  MT 
0.6889E+00 
0. 6279E+00 
0.5 59 8 E+00 
0.5065E+00 
0. 5188E+00 
0.5979E+00 
0.7140E+00 
0 .8 152E+00 
0. 8853E+00 


MAG  MP 
0.7014E+00 
0. 6230E+C0 
0 . 5393E+C0 
0.4449 E+CO 
0.4224E+C0 
C.4969E+00 
0. 6457E+  CO 
0.8003 E+00 
0.  8972E+00 


SIGTHETA 
0.1  53  5E  +  00 
C.1488E+00 
0. 1 35  5E+00 
0.1 153E+00 
0.9  112E-01 
0.6  592E-0  1 
0.4260E-01 
0.2 343 E-01 


SIGPHI 
0.1535E+00 
0 . 1523E+00 
0.1489E+00 
0. 1433E  +  CO 
0. 136  CE+OO 
0. 1272E+00 
0.1176E+C0 
0. 1074E+C0 
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C.9801E-02 
0.2164E-02 
0.1 589E-03 
0.2810E-02 
0.8790E-02 
0.1666E-01 
C.2  504E-01 
0.3278E-01 
0.3  895E-01 
0.4290E-01 
0.4426E-01 


0.9722E-01 
0.8741E-01 
0. 7831E-C1 
0.7015E-01 
0.6308E-01 
0.57175-01 
0. 5241E-C1 
0.48 795-01 
0.4625E-01 
0.4475E-C1 
0.4426E-01 
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